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SUMMARY 

A comprehensive aerodynamic analysis program based on linearized 
potential theory is described. The solution treats thickness and attitude 
problems at subsonic and supersonic speeds. Three dimensional configura- 
tions with or without jet flaps having multiple non-planar surfaces of 
arbitrary planform and open or closed slender bodies of non- circular 
contour may be analyzed. Longitudinal and lateral-directional static 
and rotary derivative solutions may be generated. 

The analysis has been implemented on a time sharing system in 
conjunction with an input tablet digitizer and an interactive graphics 
input/output display and editing terminal to maximize its responsiveness 
to the preliminary analysis problem. Nominal case computation time of 45 
CPU seconds on the CDC 175 for a 200 panel simulation indicates the program 
provides an efficient analysis for systematically performing various 
aerodynamic configuration tradeoff and evaluation studies. 
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INTRODUCTION 


This manual describes the commands and procedures for operating the 
Aerodynamic Preliminary Analysis System. It has been compiled in 
alphabetical order to simplify useage. 

Once in the system, the analyst is in a self contained interactive 
environment wherein he can input vehicle definitions, perform aerodynamic 
analysis, and display the computed results without exiting the system. It 
is however not necessary to finish an entire job at one setting. Geometry 
and analysis can be left partially complete and resumed at a later session. 
This mode has the further benefit of protecting against large information 
losses in the event of a machine failure. 

Definition of the various program functions have been made brief with 
user options cledrly" described. .The briefness is by design as it 
has been found that practice is the best teacher of the system. 

As a result of the program flexibility and freedom of choice no 
attempt has been made to define a recommended order of application. Guides 
are suggested. The analyst is urged to develop his own preferred useage 
in the process of learning the system. 

The complete program is available from COSMIC (Computer Software 
Management and Information Center, 112 Barrow Hall, University of Georgia, 
Athens, Ga 30602) as Aerodynamic Preliminary Analysis System, LAR 12404. 
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SYSTEM HARDWARE REQUIREMENTS 

Tins system, as it is available from COSMIC, was designed for operation 
on a CDC CYBER 175 and other CDC computing hardware. When compiled and 
loaded on the CYBER 175 using the included segmentation schemes, the programs 
will occupy a core space not exceeding 70k octal. 

The interactive programming has been set up to operate using Tektronix 
graphics hardware as the primary interface with the system. The Tektronix 
4014 graphics' terminal is used for keyboard input and graphics display, and 
the Tektronix 4954 digitizing tablet is used for inputing geometry 
components . Transmission speeds of less than 120 characters (1200 BAUD) 
per second reduce the graphical efficiency of the interactive programming 
and inhibit the efficiency of the system. Interested parties should 
verify the existence of such hardware and speeds at their facilities prior 
to acquisition of the program. 

Strapping options on the Tektronix- 4014 circuit boards will vary from 
installation to installation depending on computing equipment interfaces * 

The options used in the initial installation are shown on Table I for 

the TC-1, TC- 2 and Tablet interface control boards. It is recommended that 
the strapping options not be altered until after the program has been run 
once or twice since the hardware is probably correct for that 
installation. For further information on Tektronix hardware and setup, see 
Tektronix publications entitled ”4014 and 4014-1 Conprter Display Terminal 
Users”, and ”4953/4954 Graphics Tablet”. 
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TABLE I TEKTRONIX 4Q14 STRAPPING 



LF LF-CR 


TC-2 



IKfci 
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CONFIGURATION GEOMETRIC SIMULATION 


The geometry description selected for this system is the method of 
stacked sections where each section of a component contains the same number 
of defining points as indicated. 



n 


Each section can be broken into segments to improve point to point 
dependencies, reflect segments which are wetted (exposed) or un-wetted, and 
improve curve fitting on individual sections as illustrated below. The total 
number of points in any one component is limited to 500. 



SEGMENT 5 TYPE 3 and 4 

PTS 23,26 

TYPE 1 and 2 
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The system contains eight types of components for geometric construct- 
ion and analysis. 


TYPE 1: CENTERLINE BODY. Maximum number of segments per cross section 

is ten. The alignment of each section is perpendicular to the X-axis 
at construction although this can be changed after input. Components 
of this type are fuselages and centerline nacelles. 


TYPE 2: OFF-SET BODY. Has the same characteristics as a TYPE 1. This 

component is used to simulate engine nacelles , fuel pods , etc . 


TYPE 3: HALF SURFACE. Geometric quantities are based on a non-ref lected 

trapezoidal representation. 



This option is used to simulate vertical tails, end plates, and other 
components which do not require a reflected image. When a TYPE 3 surface is 
edited the view displayed is an unwrapped pictorial based on a line 
described by the points comprising the leading edge of each section. When 
a component section is listed it is displayed as a full faced cross section 
cut. TYPE 3 components (.arid TYPE 4) normally have two segments per cross 
section, an upper surface and a lower surface, but three are permitted 
to allow for airfoils with blunt trailing edges. 

TYPE 4: SYMMETRIC SURFACE. Geometric quantities are based on a reflected 

trapezoidal representation . 
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Wings are typical TYPE 4 components . The symmetric surfaces in the local 
file are scanned to select the reference first suggested for the total 
configuration. TYPE 4 surfaces are edited by displaying the component 
as it appears in the X~Y plane. 

TYPE 5: INTERFERENCE SHELLS or flat plate surface constructions. Input 

using the INTERFERENCE function or card data. 

Multiple interference shells associated with a particular slender 
body are combined into a single component . TYPE 5 constructions will 
default to TYPE 3 if a slender body (TYPE 7 or 8) is not present. 

In the EDIT function, interference shells are unwrapped and displayed 
so the overall paneling distribution can be viewed. The EDIT function 
is also used on chordwise paneling of interference shells. The LIST 
function is used to add segments. 
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TYPE 6: JET FLAPS. Input using the INTERFERENCE function. A jet flap 

is constructed in the presence of a surface component to which it is 
attached. The surface trailing edge panel chord length establishes the 
first panel of each jet flap section. Axial geometric panel spacing is 
used in EDIT and INTERFERENCE. 

Under EDIT a TYPE 6 is treated identical to a TYPE 5. 

TYPE 7: CENTER-LINE SLENDER BODY. Used in solving the isolated body 
solution, see TYPE 8 description below. 

TYPE 8: OFF-SET SLENDER BODY. Slender bodies are generated by combining 
TYPE 1 or TYPE 2 components to build up vehicle components . Bodies 
constructed of l's and 2's are designated TYPE 7. Bodies made up of all 
TYPE 2’s are designated TYPE 8. The user is limited to two slender 
bodies when performing the isolated body solution. 

Slender bodies are designed to conform to a geometric simulation 
of the input bodies which can be utilized by the isolated body program. 
The cross sectional description of TYPE 7 and TYPE 8 are 



described from bottom to top counter-clockwise with points evenly space 
about the cross section. TYPE 7 sections have 20 points , TYPE 8 sections 
have 40 points (maximum) . 
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COMPONENT NUMBERING 

The procedure for cataloging components using a numbering system is 
subject to the following rules. 

1. Conponent numbers must be greater than 50 in order to make the 
distinction between record location and component numbers when 
using the system. 

2. The viscous drag program will search for consecutive component 
numbers when a type 1 or 2 component is analyzed in order to treat 
the combined components as a single body. U£ to ten consecutive 
bodies are allowed. Bodies always start at a decimal value (110, 

300, 420, etc) while bodies to be grouped with another body have 

a unit increment from the previous component . For example components 

100, 101, 108, 109 would all be grouped together and analyzed 

as one body by the viscous drag program. 

3. Interference shells must be numbered within 10.0 of the slender body 
they belong to. Multiple- interference shells (more then one shell 
for a slender body) must sit in the working (local) file with the 
most forward shell first and the remaining shells in order behind 
it, see below. 


FILE: 

WORK 



REG 

COMP NO. 

NAME 

TYPE 

1 

700.00 

SLENDER BODY 

7 

2 

701.00 

FOR. INTE SHELL 

5 

3 

702.00 

MID INTE SHELL 

5 

L 4 

703.00 

AFT INTE SHELL 

5 


4. Type 4 components when numbered within 10.0 and stored in consecutive 
records in the work file will be combined together as one component. 
The inboard component should be followed by successive outboard 
components, see below. 


FILE: 

WORK 



.4EC 

COMP NO. 

NAME 

TYPE 

11 

400.0 

INBD WING 

— - — 7 

12 

401.0 

CNTR WING 

4 

13 

402.0 

OTBD WING 

4 
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FILE ORGANIZATION 


The geometry is stored in the permanent file using the CDC MS system 
to have direct access to any component in the file. The first record 
contains the number of components, title and the component directory. 

Records 2 through 51 are for component storage, one component per record, 
a fixed 1973 words in each record. The storage locations of each record 
and a description of each variable are given in table II. 

The variables are placed in order of importance so quantities needed 
more regularly can be accessed without reading an entire record. Each 
component is sufficiently complete that a minimum of calculation is necessary 
to display, change or analyze the geometry. 
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TABLE II COMPONENT FILE COMPOSITION 


RECORD 1 is the contents file 

LOCATION 

VARIABLE 

DESCRIPTION 

1 

NR 

NUMBER OF COMPONENTS IN FILE 

2-19 

TITLE (l)-TITLE (18) 

TITLE OF FILE 

20 

MP 

FILE UNITS FLAG 

21-70 

CCMP(1)-G0MP(50) 

SEQUENTIAL LISTING OF CCMPONENTS 
IN FILE 

RECORDS 2-51 are for component storage 

LOCATION 

VARIABLE 

DESCRIPTION 

1 

COMPN 

COMPONENT NUMBER 

2 

ITYPE 

COMPONENT TYPE 

3-6 

ANAM(l) -ANAM(4) 

COMPONENT NAME 

7 

NH 

NUMBER OF CROSS SECTION 

8 

NPH 

NUMBER OF POINTS PER CROSS SECTION 

9 

NX 

NUMBER OF SEGMENTS PER CROSS SECTION 

10 

NT 

NUMBER OF LONGITUDINAL DISPLAY LINES 

11 

IEDD 

WING DATA STORAGE FLAG 

12-21 

NMAX(1)-NMAX(10) 

NUMBER OF POINTS PER SEGMENT 

23-31 

N0Lrr(l)-N0llT(10) 

WETTED SURFACE FLAG FOR EACH SEGMENT 

32 

XMN 

MINIMUM X OF COMPONENT 

33 

XMX 

MAXIMUM X 

34 

YMN 

MINIMUM Y 

35 

YMX 

MAXIMUM Y 

36 

ZMN 

MINIMUM Z 

37 

ZMX 

MAXIMUM Z 

LOCATION 38-51 

are used for type 3 and 4 only 

38 

SW 

TRAPEZOIDAL AREA 

39 

AR 

ASPECT RATIO 

40 

TAPER 

TAPER RATIO 

41 

SWE 

SWEEP (LEADING EDGE) 

42 

DIH 

DIHEDRAL 

43 

SPA 

SPAN 

44 

CR 

ROOT CHORD 

45 

CT 

TIP CHORD 
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LOCATION 

VARIABLE 

DESCRIPTION 

46 

CB 

MEAN AERODYNAMIC CHORD c 

47 

XFCB 

SPAN DISTANCE FROM n = 0 to c 

48 

AN 

INCIDENCE 

49 

XBAR 

X-LOCATION OF l/4c RELATIVE TO 
CONFIGURATION 

50 

ME 

UNIT SYSTEM FUG 

51 

INFF 

WING DATA STORAGE FUG 

52-551 

X(1)-X(SOO) 

X-VALUES OF POINTS IN COMPONENT 

552-1051 

X(l) -Y(500) 

Y-VALUES 

1052-1551 

Z(l) -2(500) 

Z -VALUES 

1552-1571 

IXX(l) -IXX(20) 

POINTS OF LONGITUDINAL DISPLAY 
LINES 

LOCATIONS 1552-1973 are for type 3 to type 6 only 

1572 

IFLP 

NUMBER OF PANELS IN COMPONENT 

1573-1972 

IFLA(l) -IFLA(400) 

PANEL BOUNDARY CONDITION ARRAY 

1973 

INFLAP 

NUMBER OF CONTROL SURFACES IN 
COMPONENT 
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SYSTEM CURVE FITTING ROUTINES 


In order to interpolate, smooth and curve fit digitized 
sections and accomplish the other curve fitting and smoothing computations 
located throughout the system, a set of specifically designed routines have 
been set up. They are (1) linear interpolation, (2) a single cubic spline 
fit, and (3) least squares smoothing. 


The linear interpolation routine interpolates straight line segments 
evenly spaced along the set of input points. Its principal uses in the 
system include finding span references for locating sections at specified 
eta,s (2y/b), making cross section cuts in the three-view display and 
interpolating points around the outline section cuts in the SLENDER function. 


The cubic spline routine is an adaptation of a simple cubic spline which 
is point dependent at the endpoints , not slope dependent, giving the user 
spline control by the nature of his input points. This routine will also 
calculate derivatives, and do step-wise Simpson integration up to specified 
point locations. This is the principal routine used for cross section 
insertion and interpolation and segment interpolation. The interpolation is 
basically done as a function of X for TYPES 1, 2, 7 and -8 components and as 
a function of T\ for TYPES 3 through 6. The cross section cutting is done in 
conjunction with a break routine which is discussed below. 

The least- squares smoothing routine is a fixed endpoint routine 
allowing smoothing orders up to 9 but never exceeding the number of free 
points on a line segment. This routine has a built in leading edge radius 
specification for interpolating ' airfoil sections. The routine will accept 
a leading edge radius from zero or above or calculate its own radius if 
less than zero. Fuselage longitudinal lines can be smoothed with this routine 
in EDIT where they are passed through the break routine prior to smoothing. 

The break routine is designed to interpolate along specific lines on 
components breaking interpolation and starting again if a break in the 
line is detected. The break is searched for by checking the states of the 
curve behind a point relative to the changes which are about to occur at 
the next point. The routine actually begins performing between the third 
and fourth component cross sections, if there are less then four sections, 
the routine is bypassed. 
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Breaks are found in a defined line using the following algorithm i 



L-z 

Curve Illustrating Typical Break Point 
Solved for using Break Algorithm 

J 3 1 + 1 
k- 1+2 
X* i +3 

dXl » 

dyi 3 Yj -Yt 

r t = s/dx* + dyf 

Sj * /r l C sine of on^le. relative to L and i ) 

^k 5 dyj l c j C Sine of anqle relative to j and k ) 

dS k = (S k - Sj *) / dXj 

dS' k . CdS k -dSO/dS K 
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where dS’k is a check on the rate of change of the sine before point k 
relative to the rate of change in the sine about the point k. The empirical 
bias which has proved most successful to locating break points in a line 
is 


dS' k >2.45 
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FUNCTION DESCRIPTIONS 


The basic format of the system operation is to request a function 
available in the system, do what ever work is necessary in the function, 
then return to the **0K** mode of the system. For users familiar with time 
sharing system, the **0K** mode is similar to the READY or COMMAND modes 
in the two most widely used systems (TSO on IBM and INTERCOM on CDC 
respectively) . For CDC users some of the functions (CATALOG and ATTACH) 
have similar definitions to their counterparts on INTERCOM and thinking 
in terms of components in storage being the equivalent to permanent files 
on INTERCOM may ease you into the system faster. IBM users will have no 
such conflicts since similar functions yield similar results. 

The function descriptions which follow and the function dictionary 
are listed in alphabetical order. Sub- commands in the dictionary for any 
functions are listed in alphabetical order following the function to which 
they belong. Since an alphabetical description does not give the user any 
feel for overall program operation, figure 1 is provided to catagorize 
functions in terms of system usage. Also included in Appendix III are 
sample sessions illustrating cases of geometry input, component manipulation, 
analysis set-ups and results display . Using this information as a guide, 
the user may begin to operate the system. 
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ANALYSIS: Analysis provides local execution of viscous and wave drag 

solutions and set-up for flagging and processing the panel and isolated 
body solutions which are run non- interactively outside the control of the 
program operating system. The system is designed to run, and functions 
correctly, when the same Mach numbers are used for each requested 
solution. Obtaining solutions for a given configuration therefore requires 
some forethought as to the overall Mach number range to be used and then 
requesting solutions for Mach numbers covering the entire range. 

Certain solutions must be determined before others can be carried out . 
There are, consequently, some guide lines which the User should follow 
when using the ANALYSIS function. 


1) ANALYSIS is performed on all components in the local file and only 
on components in this file. All components necessary to produce the 
required solutions must therefore be attached to the local file before 
execution is initiated. 

2) The user will be informed when entering the analysis mode whether 
or not any previous solutions have been stored in the output file and what 
Mach numbers were used. If the user selects new Mach numbers all previous 
solutions are voided. Previously stored solutions are not affected if old 
Mach numbers are used. 

3) The solutions and/or Mach numbers desired on this pass are input 
by the user. Once he has set up all the solutions he enters CALCULATE 

and the computations requested are setup or executed in the proper sequence. 

4) If the user has requested 1001 suction drag he must have previously 
calculated a vortex panel solution or the wave drag due to angle of attack 
will default to zero and only wave drag due to volume will be analyzed. 

ATTACH ; Copies components from the permanent file into the local file. 

CARD: Accesses routines designed to input data in card image format. 

APPENDIX I defines the data which must be supplied to input components 
using this method. 
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CATALOG: CATALOG takes the component presently in core and places it in 

the next available record in the file specified. PERMANENT, LOCAL, 
or BOTH, are allowable parameters; PERMANENT being the default. 

COPY - The copy function is designed for copying the geometric data of a 
component into another record in the same file. This function is 
particularly useful for planar components which have been digitized with a 
large number of points for accurate definition but which must also be used 
for lifting surface solutions. The maximum number of chordwise panels 
allowed by the lifting program setup is twenty and by copying « a component and 
repaneling it with the EDIT function a well defined component can be 
brought into bounds without sacrificing the original definition. The COPY 
function can also be used to duplicate a component for use in a different 
location. This can be clumsy and the translate and duplicate sub-commands 
found in the THREE-VIEW DISPLAY function are recommended for this procedure. 

DELETE: The DELETE function is used to eliminate components from the 

permanent or local file. The file name must be provided which avoids the 
possibility that components in the wrong file will be deleted by mistake. 

If the entire local file is to be released, the PURGE function is recommended 
as the more efficient method. 


DIGITIZE: The digitize function takes components layed out as cross -sections 

on paper and input from a graphics tablet using an interactive pen or 
digitizing puck. Surface and body components are digitized and curve fitted 
or smoothed using techniques which have several important differences 
discussed below. There are also two types of components which are constructed 
using functions designed specifically for their generation, isolated or slender 
bodies and jet flaps. The digitize function also allows the user to add 
sections to previsouly digitized or otherwise input components which allows 
for replacement of simple sections with more complex ones or insertion of 
sections to improve definition . 

When digitization is completed the user is interrogated as to whether 
or not he wishes to use the spline fit or smoothing options. If the spline 
option is chosen the work is done on the component and the system returns 
to the OK mode. If the smooth option is chosen (by choosing the order of 
fit used) both the spline and smooth fits for each segment or each section 
is displayed for comparison. The order of smoothing can be changed and 
displayed anew. The user chooses at each segment which solution he wants. 
Segments containing two points are bypassed by the system. The number of 
points used in a segment is determined initially by the rniaximum number of 
points needed to define that. segment out of all sections input. This is 
fixed by first digitization and re-digitizing at a later time will not change, 
this. This rule varies for surface components and is discussed below. 
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Surface components allow (1) specific sections to be digitized or (2) 
selection of a cataloged airfoil which will be varied in thickness only. 
The section default is a 64 series airfoil. 

When the component is curve fit or smoothed after digitization a series 
of events takes place: 

1) If the smooth option is selected the user must select whether or 
not he wishes to use the leading edge radius calculated for each section or 
use a sharp leading edge. 

2) The dihedral of each section is determined based on the partial 
location of the sections in relation to one another such that the sections 
are perpendicular to a line found by the points at the leading edge of each 
section. The exception is a section of a TYPE 4 component on the center 
line which automatically receives zero dihedral. 

3) If the smoothing option is chosen the first and last Cupper and 
lower surface) segment points are summed and an equivalent number of points 
defined as : 


NP (.NFIRST + NLKST*) 

is assigned to the first and last segment. If the spline fit is chosen the 
number of points assigned is defined as 

NP = — (nfirst -*■ NUNStO 

D 

4) The scheme for curve fitting is passed on a half -cosine spread 
of the x-values spaced on the curve S along the points defining the segment 
under study. The z -values are then interplated directly against these 
x-values which are then rotated into y and z coordinates by the dihedral 
assigned to the section. This interpolation scheme provides for maximum 
definition at the leading edge of the sections where the largest slope 
variations usually take place. 

Non-planar components take the segments of each section as they are 
input and does its fitting on points evenly spaced along the curve defined 
by all the input points in that segment. No additional points are added 
to the maximum input for the segment under study. 
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DISPLAY : The DISPLAY function allows the user to either display 

orthographic projections of component combinations or three-views. Each 
option has its own set of sub -commands. 

1) Orthographic Projection 

When a single component is being displayed the user can enhance 
or retract the display lines which connect sections and produce the solidity 
in the views displayed. 

Component (s) can be displayed at different viewing angles to facilitate 
in the verification of input configurations. 

2) Three -Views 

When using the three-view option the analyst can (1) display the 
three-view, (2) cut cross-sections by selecting that function and locating 
the cross-hairs for the x- location and vertical position of the cross 
section, (3) translate components to a new location and save or duplicate 
them at the new location. When cutting cross-sections it is best to work 
from left to right to allow placement routines to work properly. 

EDIT: EDIT is a visual editing (LIST is a numerical edit function) 

function using specifically designed views of each component to insert 
sections by interpolation, delete sections, move sections up or down 
(body components) in the viewing plane, or scale sections up or down 
based on the selected chord length (surface components) . The various options 
are executed using the cross-hairs with keys on the console controlling 
the function selection. 

In addition, the surface components have a special option for chordwise 
paneling and control surface set-up. Automatic chordwise paneling for 
even spacing, half-cosine spacing, or full cosine spacing is provided 
Geometric spacing is provided for jet flaps using the panel directly in 
front of the jet flap sections to develop the spacing. The user can also 
produce his own spacing by using the cross-hairs to construct paneling. 

He can also specify control surfaces with selected keys from the console. 


EXIT: EXIT closes the perm file and plot file and ends the program. 
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FORCES : The function FORCES is used to display total forces and drag 
solutions solved for in ANALYSIS. Force data, wave drag, viscous drag, and 
drag-due-to-lift factor Cdi/Cl 2 is displayed versus Mach number. Induced 
drag is displayed as a function of lift coefficient along with any chosen 
suction variation. The characteristics of the individual displays are 
described below. The function symbol for FORCES is FORC**. 


Force data displays are requested by asking for the force and the 
parametric variable (a, 3, p, q, r, etc.) desired. Three displays can be 
shown on a single page. The coefficient combinations suggested under the 
FORCE sub -commands are the coefficients which are most often requested. 
However any combination of variables in a solution can be displayed. 

Drag due to lift is displayed versus lift coefficient and as Cdi/Cl 2 
when displayed against Mach number. Drag due to lift is presented for 
either 0% suction, 100% suction or a semi-empirical drag calculation between 
the two limits. Suction variations are displayed along with limit polars 
when this latter option is requested. 


Drag polars are calculated with suction variations based on the equation 

c * = c <w/. + > t c °«%- c D.oov.) J s=HCl) 


where S is the fractional variation of suction. The equation is applied 
at a constant lift coefficient. 


As explained in volume I, the zero percent suction drag is found from 
the integration of the surface pressures. The 100% suction drag includes the 
vortex drag from the panel solution plus the incremental far field wave drag 
due to angle of attack. In addition, an estimation of transonic divergence 
drag is added to both 0% and 100% drag polars based on the following equations, 
where Mpp is the input drag divergence Mach number: 


Mod — AAoo — 0.05 Gl. 

AAA = AA' - Aa'dc, 

= | 0.30 ( AAA+0.02) 2 4 AAA < 0 
0 I 2.50 AAA 2 + 0.00 12 , AAA 2: 0 

Also included in the 0% and 100% drag polars is any drag due to the 
basic load, i.e. wing camber, body camber, initial surface deflection. 
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Wave drag (Cqm) (far-field supersonic drag due to volume) is displayed 
versus Mach number and includes an empirical estimation for zero lift 
divergence drag in the transonic regime. The equations used to estimate 
divergence drag were discussed above in the drag due to lift section The 
laser is requested to read the FORGES section for a full explanation of the 
data reduction process. 

Viscous drag (Cpp) is displayed versus Mach number at constant pressure 
and temperature. When viscous drag is displayed, all the solutions that 
were solved are presented. 

INTERFERENCE : The INTERFERENCE function is used to build components which 

will be paneled for use in the UDP solution. These components Include 
interference shells, jet flaps, and surface constructions used for connect- 
ing components together, simulated surface fuselages for yaw or pitch 
results, etc. 

The isolated body interference shell, when declared, has it’s leading 
and trailing edge fixed by the first chord length input. Additional chords 
need only the y and z location and/or connecting components defined. 

The jet flap needs only the leading edge and connecting component declared. 

The length of each jet flap chord will be requested as a function of the 
connecting chord length. Standard panel build-up requires the leading 
and trailing edge to define the chord. Also note that when a slender body 
is not detected in the local file at UDP set-up, any interference shells 
will default to type 3 bodies. If the component used to build a jet flap 
is not present at UDP set-up a level error will result informing the user 
that the jet flap does not connect to anything and execution will continue. 

Isolated bodies and planar bodies are constructed with only one panel 
and must be edited to increase the panels . Jet flaps have two panels when 
constructed, the first set up specifically for geometric panel spacing. 

LIST : LIST is a special edit function designed for single section 

manipulation, displaying, and coordinate listing. Sections can be deleted, 
duplicated and inserted (spline-fit) however they must stay within 
the bounds of the body being, listed. Increments can be added, values changed, 
and sections of non-planar components can be rotated to simulate inlets and 
other shapes. The use of this routine requires practice, its importance 
in the system will then begin to become apparent. 
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The function symbol for LIST is 


LIST** 


and is displayed after each completed command. 

LOADS: The LOADS function is used to display chordwise and spanwise loads 

calculated with the UDP solution. It will also display the Z/c due to 
camber or thickness of con^onents set up for UDP; see figures 2 through 
4 for examples of LOADS displays. 

The function symbol for LOADS is 


PSPAN** 


PURGE : PUP-GE is used to clear the local file for executing a new case in 

ANALYSIS. PURGE sets the record value to zero and is more efficient 
then deleting all components . 

RENAME: RENAME is used to change the name of a component 


^MEMBER : RENUMBER is used to change the number of a component 


■. SAVE ^ s u f? d on a component which has been edited and updated 
to replace the old copy with new edited copy. ^ 


SLENDER: SLENDER is used to construct slender bodies from previously input 

TYPE 1 and 2 components . The components can be stacked one behind the 
other or next to each other. When components are set side by side, an inlet 
under a fuselage for example, the two parameters input under SLENDER 
(FACTOR, NPASS) have significance since they determine the bias for section 
changes; see figure S. 
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Figure 5 Cross Section Cut Before and After 
Outline Generation 

The outline routine starts at the top most point of a cross section [(point 
1 component 1, (or 1,1)]. When the outline routine gets to point . 5 of 
component 1 (1,5) it checks the distance from point (1,5) to (1,6) and 
(1,5) to (2,1) and forms the relationship: 

FACTOR. ( R[ (UO])> R C2*D] 


where R [(1,5)— ► (1,6)] is the distance from (1,5) to (1,6). If the 
relationship is true the outline skips to the second component and follows 
it as the outline until the next skip takes place (in this case from (2,11) 
to (1,10) and so on until the outline returns to the centerline or the 
first point. NPASS is the number of points passed up in the outline check 
on the section skipped to avoid back- lash problems and speed-up a lengthy 
calculation procedure. In figure 5, NPASS = 2, FACTOR =0.7 would be 
sufficient for the program to follow the section outline. 

The SLENDER function checks for inlets and exhaust nozzles prior to 
body termination and combines up to six components into a single construction 
for isolated body solutions. The user should keep in mind that each 
component input for later slender body conversion should be input in the 
same sense, i.e. either all chockwise or all counter clockwise. 

TERMINAL : The TERMINAL function allows the user to input simple geometries 

by -typing in basic descriptive parameters into the system. 
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For surface components (TYPES 3 and 4) the inputs are aspect ratio, 
area, taper ratio, sweep and dihedral, The analyst then selects the type of 
wing section to use (five different types, see figure 6) and the root 
thickness and tip thickness. 

Body components (TYPES 1 and 2) are input by selecting strategic points 
on the body (pilot station, engine face and exhaust for example) and specifying 
the x- station, cross sectional area (in units of x) and the ratio of the 
width to height of the sections. The routine performs an Emlord isolated 
body optimization on the cross sections input and interpolates one Emlord 
section in between each input section. Components constructed in this way 
have one segment per cross section and the sections input are all geometrically 
similar for the whole component. Cross sections include ellipses, half 
ellipses, rectangles and triangles. 

TIME : The TIME function keeps track of the CPU and the approximate cost 

accumulation during system execution since TIME was last called. The 
cost algorithm should be updated periodically to reflect most recent costs 
per CPU second. 

TITLE : The TITLE function perscribes a title to the permanent or local file 

which is displayed when editing and displaying components . The local file 
title is transmitted into analysis and into output files and displays. 

UNITS : The UNITS function is used to set the units of the geometry file 

and analysis. The options available are English, Meters and Centimeters, 
the default being English. English units store reference areas in square 
feet, meters and centimeters have reference areas stored in their own units. 
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COMMAND EXECUTION DICTIONARY 

The Preliminary Aerodynamic Analysis System operates in an interactive 
time sharing mode. The system will respond quickly and correctly if the 
proper input is provided. If incorrect instructions are given, the system 
will make every effort to redirect the user without falling out of execution. 

The following comments were made by experienced users and should act 
as guidelines to the help section of the manual which follows. 

1. rise of the system is the best teacher 

2. CATALOG new components after construction to avoid destroying 
the copy in core. 

3. EDIT new surface components to insure paneling does not exceed 
chordwise panel limits of 20. Since the user is limited to 200 
panels in the UDP solution dense geometric description 

should be "thinned". 
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-ANALYSIS- 


FUNCTION ANALYSIS 



OPERANDS: NONE 


ABBREVIATION: ANALY 


FUNCTION: BRINGS IN ANALYSIS FUNCTIONS. 
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-VI SCOUS - LIFTING -WAVEDRAG .CALCULATE -MACH- END- 


SUB- COMMANDS : VISCOUS, LIFTING, WAVEDRAG, CALCULATE , MACH, MDD, END 


COMMAND ^MACH NO.l, ,MACH NO. (UNDER MACH ONLY)] [ ,P (WAVEDRAG ONLY)] 


OPERANDS: VISCOUS: VISCOUS DRAG SOLUTION FLAG 

LIFTING: CONSTANT PRESSURE PANEL FLAG SOLUTION (SET-UP ONLY) 

WAVEDRAG: WAVEDRAG SOLUTION FLAG 

CALCULATE: RUN THROUGH SOLUTIONS AND/OR SETUP OF FLAGGED PROGRAM 

MACH: MACH NUMBERS TO PERFORM SOLUTION AT 

P: PRESSURE DRAG TO BE INCLUDED IN WAVEDRAG SOLUTION 


MDD: DRAG DIVERGENCE [DEFAULT = 0.90] 

END END 

ABBREVIATIONS: 


VISCOUS: V 
LIFTING: L 
WAVEDRAG: W 
CALCULATE: C 
MACH: M 
P: P 
END: E 


FUNCTION: TO SET UP AND EXECUTE PROGRAMS FOR CONFIGURATION ANALYSIS. THE 

CONFIGURATION IN THE LOCAL FILE WILL BE USED FOR ALL CALCULATIONS. 
THE LIFTING SURFACE AND ISOLATED BODY SOLUTIONS ARE DONE OUTSIDE 
THE MAIN FRAME OF THE PROGRAM DUE TO SIZE LIMITATIONS. 
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VISCOUS DRAG PROGRAM OPTIONS 


SCALE (TO FULL SIZE) : "20.0" (scale from drawing data was input from) 

1 * PLOT * NX, X(l) , . . , ,X(NX) 

2 * PLOT * NX, XB, XE 

3 * NO PLOT * NX = SI (ALL SURFACES FLAGGED) : "3" 


NOTE: 


If It 


signifies user response. Typical numbers are indicated. 
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TO EXIT OPTIONS 1 and 2 enter J0,j9,0 for answer to NX, etc. 

OPTION 1 * PLOT NX,X(1) . . ,X(NX) , (Used to remove non-wetted segments from cross sections) 


Ui 

-p* 



The cursor ( — i ) appears on screen after cross 

section is plotted. I Position cursor at one end of non- 
wetted surface line; then key in one of the record numbers 
of the components involved in the intersection. When cursor 
reappears, locate other end of non-wetted segment and enter 
record number of second component in the intersection. If 
second component already has non-wetted segment (indicated 
by dashed line) enter 0. When finished or if no segments 
need to be flagged enter zero . Each pair of entrances of the 
cursor to remove segments constitutes the removal of one 
segment from the wetted area calculation of the indicated 
components . Record numbers are those corresponding to the 
listing in the local file. Recommended NX £ 5 since the 
screen becomes too cluttered with cross sections to work 
efficiently. 


OPTION 2 * PLOT NX, XB,XE (Also used to remove non-wetted segments) 
XB * Beginning X value 
XE = Ending X value 
Procede as shown in option 1 
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Option 3 * NO PLOT * NX = 51 (ALL SURFACES FLAGGED) 


If a ll the wetted and non-wetted surfaces have been flagged by the digitizer as a 
crude approximation of the wetted area Option 3 will produce wetted area and volume 
plots directly without the arduous task of removing segments. 


After the above options are exited the program will plot the surface and volume plots of 
the configuration (see example problems in Appendix III). 


The following options are then displayed: 


* 0 EXIT * 1 SF DRAG * 2 EDIT * 

04 

on 

0 returns to or continues in ANALYSIS; 2 displays a list of x-mins, x max's for 
each component and a list of total perimeter and cross section area at each calculated 
x station, 2 then returns to Option 1 automatically. 


1 clears the page for the skin friction drag input. 
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SKIN FRICTION DRAG INPUT 

ENTER SAND GRAIN HEIGHT (KS(FT)): "0.00022" 

ENTER PRESS (PSF), TEMP (R), AND DISPLAY MACH NUMBER (0.0 FOR ALL, 999. FOR NONE) 
ENTER 0.0, 0.0, 0.0 AFTER LAST CASE (8 OR LESS) 

ENTER PRESS, TEMP, MACH (DISP) FOR CASE 1: "2116., 518., 0.6" 

ENTER PRESS, TEMP, MACH (DISP) FOR CASE 2: "0., 0., 0." 

FUSELAGE ENTER: TRANS/LENG, FLAT PLATE 1 * AXIS-BODY 2: "0.01, 2" 

WING ENTER: TRANS/LENG, CK, FLAT PLATE 1 * AIRFOIL 2: "0.01, 1.2,2" 


Oi 

ON 


NOTES: 1) Display Mach number is one for which the printed results are to be 

shown. Each Mach number from ANALYSIS input will be calculated 
for each case. 


2) KS, CK are defined in Volume I. Recommended values are given in 
Table AN-1 and AN-2 

3) Regardless of units used for input configuration inputs to the 
skin friction drag are English units. Proper conversion of 
configuration dimensions will be made by program. 

Aftet selected cases have been displayed (see Appendix III 
The options: 0 EXIT, 1 SF DRAG, 2 EDIT, will again be given 
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TABLE AN-1 


TYPE OF SURFACE 

EQUIVALENT SAND GRAIN HEIGHT KS (FEET) 

Aerodynamically Smooth 

0.0 

Polished Wood or Metal 

.00167 - .00000667 

Natural Sheet Metal 

.000013 

Smooth Matle Paint, Carefully Applied 

.0000208 

Standard Camouflage Paint, Average Application 

.0000333 

Camouflage Paint, Mass Production Spray 

.0001 

Dip -Galvanized Metal Surface 

.0005 

Natural Surface of Cost Iron 

.000833 


TABLE AN- 2 


AIRFOIL 

THICKNESS CORRECTION (CK) 

AIRFOIL t/Cj^ at 30% Chord 

2.0 

NACA 64 and 65 Series Airfoils 

1.2 







-ATTACH- 


FUNCTION ATTACH 


OPERANDS: 


ABBREVIATION: 


FUNCTION: 


ATTACH, COMP (1) ,CCMP(2) , . . .CQMP(N) 

RECORD NUMBER IN PUCE OF COMP(N) 

ALL 

ATTA 

COPIES COMPONENTS FROM PERMANENT FILE AND PUTS THEM 
IN THE LOCAL FILE. 
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-CARD- 


FUNCTION CARD 


OPERANDS: 


ABBREVIATION: 


FUNCTION: 


CARD, COMP (1) , COMP (2) , . . . .C0MP(N) 


RECORD NUMBER IN PLACE OF COMP (N) 


ALL 


NONE 


BRINGS IN COMPONENT PLACED IN THE CARD IMAGE FILE AND 
PLACES THEM IN THE LOCAL FILE. 
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-CATALOG- 


FUNCTION CATALOG 


CATALOG, [,FILeJ 


OPERANDS: 


FILE - PERM [P] [DEFAULT FILE] LOCAL [ L ] 


ABBREVIATION: CATA 


BOTH 


FUNCTION: 


CATALOG COMPONENT IN CORE INTO REQUESTED FILE 




-COPY- 


FUNCTION COPY 


OPERANDS: 


ABBREVIATION: 


FUNCTION: 


COPY, COMP, COMPN ,FILE 


COMP - COMPONENT TO BE COPIED 

COMPN - COMPONENT NUMBER TO BE COPIED TO 

FILE - PERM [p] [DEFAULT FILE] , LOCAL [l] 

NONE 


COPY A COMPONENT IN REQUESTED FILE TO A NEW RECORD 
IN THE SAME FILE. 
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-DELETE- 


FUNCTION DELETE 


DELETE, FILE, CCMPl HCQMP2 , . . . [,ALL]J 


OPERANDS: FILE- PERM [P] , LOCAL [L] (MUST BE SPECIFIED) 

C0MP1 OCMPONENT(S) TO BE DELETED 
ALL - DELETE ALL COMPONENTS 

ABBREVIATION: DELE 

FUNCTION: REMOVE UNWANTED COMPONENTS FROM FILES 
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-DIGITIZE- 


FUNCTION DIGITIZE 


DIGITIZE QcCMPN] [,AXIS] 


OPERANDS: COMPN: RE -DIGITIZE COMPONENT COMPN . DEFAULT IS NEW COMPONENT 

AXIS: NEW AXIS SYSTEM WILL BE REQUESTED. DEFAULTS TO AXIS 

ON FIRST CALL TO DIGITIZE. DEFAULTS TO PREVIOUS AXIS 
ON SUBSEQUENT CALLS 


ABBREVIATION: DIGI 


FRUNCTION: 


DIGITIZE COMPONENT FROM DRAWING OF THE DIGITIZING TABLET 




-DIGITIZE- 


SYMBOL 


® 

1 ( 
L 



CR 


’’QUANT’ 


DESCRIPTION 

Point input from tablet 
Info in bracket would not appear on screen 
Indicated infoimation is on graphics screen 
Tablet input using digitizing pen or puck 
Carriage return 

Position of point to be input from tablet 
QUANT is a user input 


44 



-DIGITIZE- 


VARIABLE 

XO 

XMAX 

TYPE 


NSEG 

NOUT(I) 


TOC 


DESCRIPTION 

: MINIMUM X VALUE FROM DRAWING 
: MAXIMUM X VALUE FROM DRAWING 

:1 - NON- PLANAR CENTERLINE (EG. FUSELAGE) 

2 - NON-PLANAR OFF-SET (EG. ENGINE POD) 

3 - HALF PLANAR (EG. VERTICAL TAIL) 

4 - FULL PLANAR (EG . WING) 

: NUMBER OF SEGMENTS/SECTION 
: WhTi'EU- UNWETTED SEGMENT FLAG 
1 - WETTED 
-1 - UNWETTED 
t/C MAX OF AIRFOIL 

TOC < 0 PROGRAM ASSUMES USER IS INPUTTING 
AIRFOIL SECTION 

TOC > 0 PROGRAM USES LAST INPUT AIRFOIL SECTION 
SCALED TO NEW CHORD AND TOC 
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DIGITIZE 




-DIGITIZE- 



JZ-k 
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617 



DIGITIZE 



TOC >0 


POINTS (a)-( Jo) Se.GrWVeNT 1 
POINT'S. (TiV(T&) se&w\£NT 2 


® © 



□ P 


S \ck 


W~k 


*4 HE XT. HOOP X 2. RCPCAT HOOP* 3 END BXY 



IST-A 


INPUT OF TYPE 3 t4 S ECTIONS 
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-DIGITIZE- 


NOTES ON DIGITIZE 


1) XO and XMAX do not have to extend the length of the tablet or 
encompass the entire configuration. They are simply two convenient 
points on the drawing. 

2) The XO in the plan view and profile view must be the same X-station. 

3) When a change in segments occurs, the user inputs two points. If two 
or more segments occur at the same point (i.e. segments of zero length 
such as would occur on the cross section at the end of a canopy. See 
below) > the user should input a number of identical points equal 

to the number of segments, i.e. if three segments meet at the same point 
three points are put there., 



The first and last points on a section do not require multiple points 
unless multiple segments converge there. 
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4) If the default airfoil (64A0XX) is used it will remain the reference 
airfoil unless changed by the input of a new airfoil section. The latest 
airfoil section always becomes the reference section. The default section 
can only be used if the planar component has two segments. 

5) There are two leading edge radius algorithms in this program, based on a 
polynomial representation of the airfoil leading edge, in 
conjunction with the least -squares smoothing routine. The first method 
calculates the leading edge radius based on the input points nearest the 
leading edge . The second method is based on the 5-digit series airfoil 
radius: 


R l.EL. - I.OI** (Vc max ) Z CHORD ({.4 -Xmax /c ) 


where *max/c is the fractional chord location of the maximum airfoil 
thickness. 


6) The multipoint routine used to input points from the graphics tablet has 
two inport ant keyboard features, S \ GR ] and R \CR]. The S [CR] tells 
the program that the cross-section is finished and the points should be 
entered for calculation. R ^CR^ tells the program that an input error, 
lost track of where I was, etc., occurred and it should start at the 
beginning of the loop again. 

7) The maximum number of input points allowed per cross section is forty. 
This number includes the reference POINTS put in for a cross -section* 
such as indicating the center line for a body section or chord markers 
for an airfoil section. 
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-DISPLAY- 


FUNCTION DISPLAY 



OPERANDS: FILE: [PERM [p] DEFAULT FILE] [LOCAL [l]] 

VIEW: [ORTHOGRAPHIC [q] DEFAULT 1 [tHREE-VIEW [t]1 

OQMP(N) : RECORD NUMBER IN PLACE OF COMP(N) 

DEFAULT IS COMPONENT IN CORE 

ANG: YAW,PITCH,ROLL. VIEWING ANGLES (ORTHOGRAPHIC ONLY) 

DEFAULT ANGLES (- 65 , 25 , - 25 ) 

ALL: ALL COMPONENTS IN FILE 


FUNCTION: DISPLAY COMPONENT (S) IN ORTHOGRAPHIC PROJECTION OR 

THREE VIEW 
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-DISPLAY- 


Orthographic Projection: 

General: Components displayed using the orthographic projection can be 

rotated and displayed from any angle. The order of angle specification is yaw, 
pitch, and then roll. Each view is scaled to fit the entire screen. 

Single Component Display: When a single component is displayed the options 

are: . 

* 1 VIEW * 2 ENHANCE * 0 END 

The view option is self explanatory, as is end. The enhance option controls 
the number and location of the longitudinal lines connecting each cross section 
which give it a solid body appearance. The enhance option produces the following 
response: 

* PRESENTLY LINES AT : 1 5 10 

* ENTER NO. ADDITIONAL LINES, LINE LOCATIONS * 

The print out informs you that all point number one T s in each cross-section have 
been joined with a display line, as are all the 5's and 10 T s. Twenty such lines 
are allowed. The user simply enters the number of lines he is going to input then 
a positive point location to add a line or a negative location to remove one. 

Single components are only displayed in half view and a little block printed 
in the upper right hand comer. It is necessary to display newly input planar 
components at least once before cataloging them to fill in any blanks in the title 
block. 

Multiple Body Displays: When more than one component is displayed (they can 
be the same component) both sides of the configuration are displayed. The only 
option available are the view and end option . 

Three ^View Displays: 

The three view display offers two options, component translation, and 
cross section display, refer to figure DI-1 the three-view schematic. 
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-DISPLAY- 


Component translation notes: 

1) Component number requested is the integer value next to the component 
name as indicated (in the schematic component 3, the VERTICAL, was selected, for 
translation) 

2) Translation modes 1, 2 and 3 are two step operations requiring a 
"from point" and a "to point ". The selected mode is entered with the "from 
point". Successive moves are possible in any of the first three modes and 
terminates with a 5 input on a "from point". 

3) Translation mode four is a one-step process, the values of dx, dy, 
and dz are requested in the upper left hand corner under the component number 
line. The same results can be achieved as in modes 1, 2 and 3 and all end with the 
inputs requested in the lower right hand box of figure DI-1. 

4) A saved component is placed in the location of that component in the file 
the initial display was made from (i.e. perm or local). A duplicate of the 
component (translated to its new location) will likewise be filed in the displayed 
file in the next available record location. A new number 0.10 higher than the old 
number is assigned to this new component. 

5) You can continue on the same display or redraw the display if it has become 
too crowded. Redraw should not be used to just display since it assumes more 
translations are to taice place. 

Display Cross-Sections: 

Cross section cuts are displayed by locating the graphics cursor at the 
desired x location and locating the y-cursor where the y = 0.j9 line is desired. 
Entering a non- zero value from the keyboard produces a section at that location. 

A zero ends the routine and pages the screen; therefore a copy should be taken 
before hitting zero. 
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START 



• TRANSLATE I ♦ CROSS SECTIONS 2 ♦ END 0 








QRVUICS CUR SR 



FIGURE DI-1. THREE VIEW INSTRUCTIONS 








-EDIT- 


FUNCTION EDIT 


OPERANDS: 


ABBREVIATION: 


FUNCTION: 


EDIT [, FILE] [,COMPN] [,COMPREF] [,ROLL ,ANG] 


RECORD NUMBER IN PLACE OF COMP 

LEAVE OUT COMPN TO EDIT COMPONENT IN CORE 

COMPREF - USED ONLY FOR OFFSET BODY. 

DEFAULT IS NONE 

ROLL, ANG - USED ONLY FOR OFFSET BODY 
COMPONENTS, WHERE ANG IS THE VALUE' OF ROLL DESIRED, 
DEFAULT VALUE OF ANG IS 0 DEGREES 

FILE - [PERM [P], DEFAULT FILE], [LOCAL [L] ] 

NONE 


VISUAL EDIT FUNCTION. COMPONENT EDITING IS DONE USING 
DISPLAYED VIEWS OF THE COMPONENT. 
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-EDIT- 


BODY COMPONENT EDITING 

Figures ED- la and ED- lb illustrate how the edit subroutine displays 
a body component. The options available are summarized in the upper right 
hand comer of the figure ED- lb and reproduced below. 


INSTRUCTIONS 

POSITION SPAN STATION 
I: INTERPOLATE STATION 

D: DROP STATION 

C: CAMBER STATION 

INITIAL Z 
FINAL Z 

Y: SMOOTH Y 

Z: SMOOTH Z 

B: SMOOTH BOTH 

SPACE: END STATION INPUT 


All options are initiated by locating the vertical graphics cursor 
(thin vertical line) at the fuselage station desired and entering the in- 
dicated key on the console. With the exception of options I and C only one 
keyboard entry per station is required. 

EDIT OPTION I: Interpolates and enters a section at the indicated fuselage 

station. Duplicate stations are not allowed and will result in multiple 
rings of the "bell” on the console. (Figure ED-lc STEP III). 

EDIT OPTION D: Section chosen will be deleted from the geometry description. 

If the station indicated does not correspond to an available section the 
multi-bell function will sound. (Figure ED-le STEP III). 

EDIT OPTION C: The vertical graphics cursor is positioned at the station 

to be cambered (it must be an available station) . The horizontal cursor 
is moved to the initial Z location and a C is entered on the keyboard. The 
graphics cursor will go off and a triangle will appear at this Z location. 
When the graphics cursor appears again the final Z is located and a G is 
again entered on the keyboard. A triangle will appear at the final Z 
location and the new section will be displayed, just to the right of the 
old section, for. comparison. The fuselage stations of the new sections 
are the same as those of the old sections (Figure ED-lc STEPS I and II). 


58 



EDIT OPTIONS Y, Z, and B: When the user selects these smoothing options, 

the system applies a piecewise smoothing algorithm of specified order. 
Smoothing will be done on Y- values if Y is selected, Z- values if Z is 
selected, or both Y and Z if B is selected. The system will request the 
maximum smoothing order at the lower left hand comer of the page. The 
system will reduce the order if the number of free points is less than the 
maximum order on any line segment. 

EDIT OPTION SPACE: Depressing the space bar ends the section editing 

routine and returns control back to the main system. 


59 



-EDIT- 



FIGURE ED-1. SIMPLE FUSELAGE BUILD-UP 
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FIGURE ED-1 (CONTINUED) 
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-EDIT- 



FIGURE ED-1 (CONCLUDED) 
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-EDIT- 


SURFACE COMPONENT EDITING 

The edit options for all surface components are the same, therefore the 
wing and vertical tail will be illustrated using the wing. The distinction 
made by edit between full surface and half surface components is the manner 
with which it handles the spanwise representation of either component. Full 
planar components are viewed in an X-Y plane or water -plane perspective. 

Half surface components are ’’unwrapped” and displayed in a true view, i.e. a 
curved half planar component would be edited spanwise along an s coordinate 
where s is perpendicular to the X axis and is measured from the root along 
the surface to the tip of the component. 

The wing input is shown in its initial display in figure ED- 2a. 

Because the number of spanwise lines which can be displayed for any component 
would confuse the display, only a few selected spanwise lines are used by the 
program to represent the shape of the component in question. This is also 
true of longitudinal lines displayed for body components . 

When the user selects the edit option the program presents him with a 
view of the present status of the wing, its shape and its paneling figure 
ED- 2b. The paneling on surface components is important if the user is 
interested in the lifting surface option of the program. The edit options 
available here are: 

* OPTIONS: NO CHANGE 0 * CHORDWISE FIT SECTIONS 1 * CHANGE SPAN STATIONS 2 * 

« 

If the user intends to adjust the span stations, they are done first and 
option 2 is selected. 

The instructions for planar edit appear in the upper right hand corner 
of the screen, figure ED- 2c, and are reproduced below 


INSTRUCTIONS 

POSITION SPAN STATION 
I: INTERPOLATE STATION 

D: DROP STATION 

E; CHORD ELONGATION 
LOCATE L.E. 

LOCATE T.E. 

SPACE; END STATION INPUT 
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-EDIT- 


The graphics cursor will appear and section editing can begin. Options I, 

D, Space Bar are the same, as for non-planar components and will not 
be repeated here. The hash marks at the bottom of the screen mark every 5% 
true span for the component as guides to inputing new sections. 

Option E is for chord elongation or shrinking, requiring two inputs 
by the user, one for the leading edge and one for the trailing edge 
of the chord. The vertical cursor is located at an available section and the 
new leading edge of the section is located with the horizontal cursor 
(the multi -bell will sound if the station chosen is not available). An E 
is entered on the keyboard, the graphics cursor goes off and a triangle is 
drawn at the input point, figure ED- 2c STEP I. The process is then 
repeated for the trailing edge and the bell will ring once when the process 
is complete, figure ED- 2c STEP II. The elongation routine scales the 
thickness in direct proportion to the old and new chord lengths. The planform 
modification process performed on the example wing is shown in figure ED- 2d 
where the new leading and trailing edge points are numbered in the order in 
which they were input. Note, however, that although this user made an orderly 
advance from wing root to tip, the order in which the sections are changed, 
dropped, or added is not important in this section of the edit routine. 

When the user ends this section of edit the progrm again offers the choice 
of chordwise paneling the geometry: 

* M3 CHANGE 0 * CHORDWISE CONSTRAINTS 1 * EVEN SPACING 2 * 

HALF COS SPACING 3 * FULL COS SPACING 4 * 

Examples of choices 2 thru 4 are shown in figure ED- 3a, b, and c. Compare 
these with the aribtrary (chordwise constraints) result in ED- 3d. The 
response, by the computer to choices 2 thru 4 is: 

*‘ENT1R NO. OF CHDWSE PANELS * 

* PRESENT NUMBER = 10 * 

? 

to which the user responds with the desired number of chordwise panels. The 
program will then plot a final view of the edited component showing the 
paneling for the component and ends the editing mode, as in figure ED-4f . 
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FIGURE ED- 2. SURFACE COMPONENT SECTION EDITING 
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FIGURE ED-2. (CONCLUDED) 
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FIGURE ED-3. EXAMPLES OF CHORDWISE PANEL DISTRIBUTION METHODS 
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-EDIT- 


CHORDWISE CONSTRAINTS 

The user may elect to input his own chordvd.se constraints by selecting 

option 1. A chordwise constraint point is a user input point defining locations 
on any given section where the user requires specific boundaries to pursue the 
desired computational results. Constraints are imposed, for example, to produce 
arbitrary chordwise paneling or to indicate control surfaces on the component . 

A constraint line always runs from root to tip along a user described path. 
Constraint sections run between two constraint lines, the leading and trailing 
edges of the component being automatic constraint lines. 

Some general rules concerning flap generation are: 

1) Flaps are defined along a constant line input at the trailing edge 
of the surface . For' example, a leading edge flap would be defined along 
with first input constraint line while a trailing edge flap would be defined 
after the trailing edge has been indicated. 

2) A flap definition extends from the constraint line being input, forward 
to the previously input constraint line. 

When choice 1 is selected, figure ED-3a 9 a shell of present surface 
definition is plotted with the instructions in the upper right hand comer, 
figure FED- 3b and reproduced below: 

INSTRUCTIONS 

C: CONSTRAINT 

I: INTERMEDIATE FLAP LINE 

F: BEGIN FLAP 

E: END FLAP 

T: TRAILING EDGE 

R: REFERENCE LINE 
2ND POINT 

/: END CONSTRAINT SECTION 

The chordwise hash marks are every 5$ for reference purposes. The user is 
prompted here for the type of spacing desired in the first constraint section, 
Step 1, figure ED- 3c : 

* CONSTRAINT NO. : 1 ** EVEN SPACING 2 * COSINE SPACING 3 **, 

0 ) 

where the cosine spacing here is half cosine spacing. The program then 
requests the number of panels to be placed in this constraint section. 

Step 2: 
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-EDIT- 


* NO. CHRDWSE PANELS THIS CONSTRAINT * 
7 



The graphics cursor now appears to start constructing the first 
constraint. Note that Figure ED- 3 illustrates the steps necessary to 
produce the work which is shown. The constraint line will end 
automatically when the maximum span station is selected or when a / is input 
from any section. When a slash is encountered the x/c and flap designation 
information will be derived from the previous input, i.e. / can never be 
the first input, for a constraint line. 

In figure ED- 3c the user was mainly interested in achieving close 
spacing at the leading edge and constructed the first constraint line by: 

1) Locating the graphics cursor at the 0.25 x/c location of the most 
inboard chord and entering a.c. 

2) Leaving the graphics cursor alone when it re -appears and entering 

a /. 


The constraint line is then drawn in solid and the computer generated 
intermediate lines are drawn with a double dashed line. 

The second constraint is constructed, figure ED- 3d by indicating 
an even spacing of 5 panels. The horizontal cursor is then located at the 
0.80 x/c point of the fourth outboard section. The constraint points to be 
input define the leading edge of the aft control surfaces on the wing. They 
are input as follows: 

1) Starting inboard at section 1 a C is entered at sections one through 
4 leaving the horizontal cursor fixed throughout. 

2) Following the C input at the fourth station a /is entered and the 
constraint is ended. (Note that outboard of section 4 constraint line 2 is 
at a constant x/c of 0.8, the same as it was at 4). 

To finish the wing with constraint 3, the user indicates cosine spacing 
with, 2 panels, figure ED-3e. Since this constraint is along the trailing 
edge a T is entered when the graphics cursor comes on, indicating the trailing 
edge, releasing the need for the user to locate anything with the horizontal 
cursor. The user now proceeds, using the vertical cursor only, to locate the 
partitioning of the trailing edge control surfaces as follows: 

1) The first section is entered with a C 

2) The second section is entered with an F indicating that a flap is 
starting. 
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-EDIT- 


3) The third section is skipped and an E is entered at the fourth 
indicating that a control surface is to be entered between sections 2 and 
4. 


4) An E is entered at section 5 and a second flap is entered by the 
program between 4 and 5 (An E following another E indicates that a new 
flap should be entered between these two span stations) . 

Similar operations are done at sections 7 and 9 to indicate two 
additional control surfaces from 5 to 7 and 7 to 9. After the E is entered 
at section 9, a slash is entered and the program fills in the remaining 
constraint section on the wing and ends this section of edit by ringing 
the bell and pausing with a PCC for a copy, if desired. The program then 
proceeds to plot a final view of the component showing all paneling. 
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FIGURE ED- 4. SURFACE COMPONENT CHORDWISE CONSTRAINT EDITING 
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-EXIT- 


FUNCTION EXIT 



OPERANDS : NONE 


ABBREVIATION: NONE 


FUNCTION: CLOSE FILES AND EXIT FROM SYSTB1 
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-FILES- 


FUNCTION FILES 


OPERANDS: 


ABBREVIATION : 


FUNCTION: 


FILES [,FILE] 


FILE - PERM [P] [DEFAULT FILE] 
LOCAL [L] 


FILE 


LIST OUT COMPONENTS STORED IN EITHER PERMANENT 
OF LOCAL FILE 
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-FORCES - 


FUNCTION FORCES 


OPERANDS: 


ABBREVIATION: 


FUNCTION: 


FORCES 


NONE 


FORC 


DISPLAY TOTAL FORCES FROM PROGRAM CALCULATIONS. AT 
LEAST ONE ANALYSIS FUNCTION MJST HAVE BEEN COMPUTED 
TO PRODUCE REASONABLE RESULTS. 
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-FORCES SUB-FUNCTIONS- 


OPERANDS: 



COM, COM2, COM3: 


CLO: 

CMD: 

CMA: 

CLA: 

CMCL: 

CLP: 

CLQ: 

CMP: 

CM J: 

CLDS1* 

CMDS1* 

CLDS2* 

CMDS2* 

CYB: 

CNB: 

CRB: 

CYP: 

CYR: 

CNP: 

CNR: 

CRP: 

CRR: 

CYDA1** 

CNDA1** 

CRBA1** 


CDP 

CDM 

K0P 

KLP 

KS# 

SA: 


SO: 


C L 0 

CMo 

CMa 

CLa 

dOn/dCL 

CL p 

CLq 

S 

Qnq 

9CL/86SX 

9Qn/ 96Si 

3.0/ 96s 2 

9Qn/96S2 


CYg 

p S 

% 

Cy| 

Cnp 

Cn r 

Cjtp 

C *T 

3 V 3 «ai 

3C n/ 3 6Ai 

3Cji/ 3 5Ai 


VISCOUS DRAG 
WAVE DRAG 

c Dj/ c L 2 0% SUCTION 

CWCl? 100% SUCTION 

c d|/Cl 2 USING SUCTION VARIATION * 


ABSCISSA SCALE FACTOR 
0.0 FREE SCALING. DEFAULT 0.0 

ORDINATE SCALE FACTOR 
0.0 FREE SCALING. DEFAULT 0.0 
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ABBREVIATION: 


NONE 


FUNCTION: FUNCTIONS SPECIFIED ABOVE ARE ALL PLOTTED VERSUS MACH 

NUMBER. THREE DERIVATIVES MAY BE SPECIFIED PER DISPLAY. 


K0P,K1P,KS#,CDP,CDM ARE PLOTTED ALONE 


*DS1 REFERS TO SYMMETRIC FLAP 1 
**DA1 REFERS TO ANTI- SYMMETRIC FLAP 1 
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-FORCES ,SUB-FUNCTIONS- 


OPERANDS: 

COM: 


MACH: 

SA: 

SO: 

FUNCTION: 



ALL ARE VRS 

a 


CDI0 

* Cdi 

01 SUCTION 

CDI100 

- CDi 

100% SUCTION 

CDI# 

- CDi 

SUCTION VARIATION 


WHERE # IS SUCTION CURVE 


SUC# - SUCTION CURVE # 

CDT0 Tnv] - CUroiAL 0* SUCTION 

cnri0j3(^NvJ- cdtotal 100 % suction 

CDT# [NV] - C0TOTAL WITH SUCTION VARIATION # 

NV - VISCOUS DRAG CASE NUMBER [DEFAULT-l] 

MACH NUMBER, MUST BE AN INTERPOLATED MACH NUMBER 

ABSCISSA SCALE FACTOR 
0.0 FREE SCALING . DEFAULT 0.10 

ORDINATE SCALE FACTOR 
0.0 FREE SCALING . DEFAULT 0.0 


DISPLAY DRAG POLARS AND SUCTION VARIATIONS VERSUS LIFT 
COEFFICIENT. 
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-INTERFERENCE- 


FUNCTION INTERFERENCE 


INTERFERENCE, [,FILE]C0M REC1 .COM2 REC2 ...COW RECn 


OPERANDS: 


COM1: COMPONENTS TO BE USED FOR REFERENCE .AT LEAST ONE 
COMPONENT MUST BE SPECIFIED 
[RECl] RECORD NUMBER MAY BE USED INSTEAD OF 
COMPONENT NUMBER 


FILE: [PERM [p] , DEFAULT] , [lOCAL [l]1 
ABBREVIATION: INTE 


FUNCTION: USING SELECTED COMPONENTS THE OPERATOR CREATES INTERFERENCE 

SHELL, JET FLAPS, AND FLAT PLATE GEOMETRIES. 
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-INTERFERENCE- 


NOTES: 

1) Panel building with this function is a recurrent two-step process 
requiring at least two passes to produce a geometric figure constituting a 
planar-body. The recurrent process takes place in steps II through IV as 
depicited in figure IN-1. 

2) The three types of boundary conditions, i.e., interference, standard, 
and jet flap, are treated somewhat differently with the graphics cursors as 
explained below. 

Interference lines: End points are determined by the values indicated by 

the vertical cursor for the first input line. There after the vertical cursor 
can be moved aside and only the horizontal cursor is used. 

Standard lines: Ends points of standard lines must be located for each line 

using the vertical cursor. 

Jet Flap lines: Only the leading edge of the line must be located and each 

line must attach to a displayed chord. The length of the line is requested as a 
multiple of the attachment chord length (Step IV for jet flaps only in figure 
IN-1). 

3) Components in the display not qualifying for planar attachment are 
dashed in and can be used as reference guides only. Component such as fuselage, 
nacelles etc. are used for this purpose. 

4) Lines are input by locating the graphics cursors at the X- leading edge, 

Y location and entering an attachment component number (component numbers to use 
are printed at the approximate centroid in a full view of each component ; for 
example wings in the top view, verticals in the side view) to input the point. 
When the cursors appear again locate the X- trailing edge, Z location and input 
the point by entering the boundary condition flag (I, E, or J). Do not mix 
boundary conditions. After the X, Z point is entered the bell will ring and a 
line will be drawn in the top and side views based on the points received. Jet 
flaps draw a line of representative length in the top view only. 

5) Interference shells must be numbered within 10.0 of the slender body 
they are associated with. 

Interference shells numbered within 10.0 will be combined into one component 
in the lifting surface solution. 
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-LIST- 


FUNCTION LIST 


OPERANDS: 


ABBREVIATION: 


FUNCTION: 


LIST [FILE] [COMP] 


RECORD NUMBER IN PLACE OF COMP 

LEAVE OFF COMP TO EDIT COMPONENT IN CORE 

FILE: [PERM [p] DEFAULT] , [LOCAL [l]] 


NONE 


SECTION LISTING AND EDITING. SECTIONS CAN BE DISPLAYED, 
INSERTED, DELETED. SPECIFIC POINTS CAN BE CHANGED AND 
SECTIONS CAN BE ROTATED (BODY COMPONENTS ONLY) . 
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-LIST / A,B,BA,DB- 


SUB- COMMAND A, 


OPERANDS: 


ABBREVIATION: 


FUNCTION A: 


B, DA, AND DB 



SEC - SECTION NUMBER. DEFAULT: CURRENT SECTION 

=VAL - VALUE OF ANGLE FOR ROTATION 
+ COUNTER-CLOCKWISE 
- CLOCKWISE 

(MUST HAVE SIGN IN FRONT OR DECIMAL 
POINT IN VAL) 

N1,N2 - START POINT AND END POINT 

IF N1 IS GREATER THAN N2 ROTATION TAKES 
PLACE ON POINTS NOT BETWEEN N1 AND N2 

N3 - POINT TO ROTATE ABOUT, DEFAULT IS N1 

V - VERIFY INPUT VALUE 

DEFAULT IF N1,N2,N3 ARE LEFT OFF: 

ALL POINTS IN SECTION ARE ROTATED ABOUT 
POINT 1. 


NONE 


INCLINE BODY SECTIONS TO SHAPE INLETS, AND NOZZLE 
CONTOURS. SECTION RETAINS SAME CROSS-SECTION IN Y-Z PLANE 

A = DA = INCLINE ABOUT Y-AXIS 

B = DB = INCLINE ABOUT Z-AXIS 
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-LIST / A,B,DA,DB- 


FUNCTION B: TWIST AIRFOIL AND ADJUST DIHEDRAL 


A = DA = TWIST ANGLE, ADDED TO WHAT EVER TWIST IS 
ALREADY PRESENT 

B = DB = ABSOLUTE DIHEDRAL ANGLE OF AIRFOIL. DB =0.0 
RETURNS SECTION TO ZERO DIHEDRAL 
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“LIST / DELETE- 


SUB- COMMAND DELETE 



OPERANDS: 


N1 - SECTION TO DELETE 

N2 - DELETE SECTIONS FROM N1 TO N2 


V - VERIFY INPUT SECTION 


ABBREVIATION: DE 


FUNCTION: 


DELETE SECTIONS FROM COMPONENT GEOMETRY. 
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-LIST / DUPLICATE - 


SUB-COWAND DUPLICATE 


DUPLICATE-VAL ,N1 [, V ] 


OPERANDS: 


=VAL - STATION. FOR NEW SECTION 
x FOR NON-PLANAR 
Y FOR TYPE- 4 (OR WING) , Y > 1.0 
r) FOR TYPE-3, 4, 5, 6 

N1 - STATION TO BE DUPLICATED 

V - VERIFY INPUT SECTION 


ABBREVIATION: DU 


FUNCTION: DUPLICATE STATION N1 AT LOCATION VAL. 
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-LIST / DX,DY,DZ- 


SUB- COMMAND DX, 


OPERANDS: 


ABBREVIATIONS: 


FUNCTION: 


DY, AND DZ 


[sec,]dx =VAL [,Nl][,N2] [,v] 


SEC - SECTION NUMBER. DEFAULT: CURRENT SECTION 

=VAL - VALUE TO BE ADDED TO X 

(MUST HAVE "=" SIGN IN FRONT OR DECIMAL 
POINT IN VAL) 

N1 - POINT TO HAVE INCREMENT ADDED 

N2 - ALL POINTS FROM N1 TO N2 ARE TO HAVE 
INCREMENT ADDED 

V - VERIFY INPUT VALUE 

DEFAULT WHEN Nl, N2 ARE LEFT OFF IS "ALL" 


NONE 


ADD INCREMENTAL VALUE TO X, Y, OR Z QF A GIVEN 
SECTION USING SUB-OMAND DX,PY, OR'DZ, 
RESPECTIVELY 
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-LIST / INSERT - 


SUB-COMMAND INSERT 


INSERT =VAL [ ,V ] 


OPERANDS: VAL - X STATION OF NON-PLANAR 

Y STATION OF TYPE- 4 (OR WING), Y > 1.0 
n STATION OF TYPE-3, 4, S, OR 6 


V - VERIFY INPUT SECTION 


ABBREVIATION: IN 


FUNCTION: CURVE FIT A NEW SECTION AT INPUT STATION. SPLINE FIT 

ROUTINE USES AN ALGORITHM DESIGNED TO FOLLOW BREAKS IN 
THE COMPONENT. 
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-LIST / LIST- 


-SUB-OM1AND 


OPERANDS: 


ABBREVIATION: 


FUNCTION: 


LIST 


[sec,] list [,rot] [,sp] [,sd] 


SEC - SECTION NUMBER. LEAVE OFF FOR CURRENT 
SECTION 

ROT - ROTATE SECTION ABOUT Z-AXIS 

(MUST HAVE DECIMAL POINT IN VALUE) 
DEFAULT IS FACE ON VIEW 

SP - SUPPRESS PRINT OUT OF POINT X, Y, AND 
Z VALUES .DEFAULT . IS PRINT 

SD - SUPFRESS DISPLAY OF SECTION . 

DEFAULT IS DISPLAY 


LIST OUT SECTION DISPLAY AND X, Y, Z PRINT-OUT. 
MULTIPLE POINTS ARE NOTED BY BREAKS IN NUMBERING 
ON THE SECTION DISPLAY. 
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-LIST / SAVE- 


SUB- FUNCTION 


OPERANDS: 


ABBREVIATION: 


FUNCTION: 


SAVE 



CQMPN 


] 


CQMPN - SAVE AS NEW NUMBER . 
DEFAULT IS COMP 


S 


SAVE WORK DONE IN FILE COMPONENT WAS TAKEN FROM. 
NEW COMPONENTS ARE SAVED IN PERMENANT FILE. 
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-LIST / X,Y,Z- 


SUB-COWANDS 


OPERANDS: 


ABBREVIATION: 


FUNCTION: 


, Y, AND Z 


[SEC]X[=VAL] [,Ni][N2] [,N3] [,V ] 


SEC - SECTION NUMBER. DEFAULT: CURRENT SECTION 

=VAL - REPLACEMENT VALUE OF X MUST HAVE 

SIGN IN FRONT OR DECIMAL POINT IN 

VAL) 

N1 - POINT TO BE GIVEN NEW VALUE 

N2 - ALL POINTS FROM N1 TO N2 ARE TO BE 
SET TO VAL 

N3 - ALL POINTS FROM N1 TO N2 ARE TO BE SET 
THE VALUE OF POINT N3 

V - VERIFY INPUT VALUE 

DEFAULT WHEN Nl, N2, N3 ARE LEFT OFF IS "ALL" 


NONE 


REPUCE PRESENT VALUES OF X, Y, OR Z OF A 
GIVEN SECTION USING SUB- COWAND X, Y, OR Z, 
RESPECTIVELY. 
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-LOADS- 


FUNCTION LOADS 


OPERANDS: 


ABBREVIATION: 


FUNCTION: 


LOADS 


NONE 


LOAD 


DISPLAY CHORDWISE AND SPANWISE DISTRIBUTIONS 


THE SET-UP AND/OR THE SOLUTION TO TOE LIFTING SURFACE 
FUNCTION MUST BE EXECUTED BEFORE THIS’ OPTION IS 
SPECIFIED. 
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-LOADS SUB FUNCTIONS- 



OPERANDS: 

CCH: CP - Cp UPPER AND LOWER VRS X/C 

CPN - Cpnet VRS X/C 

U - U/V« VRS X/C 

Z/C - Z CAMBER/C yRS x/C 

Z/T - Z THICKNESS/C VRS X/C 

CLC - C t C/CavG VRS SPAN (n) 

CN - C n VRS SPAN 

CNC - CnC/CAVG VRS SPAN 

CM - X/C C.P. VRS SPAN 

END - [E] END 

ADD: ADD ANGLES LISTED TO THOSE PRESENTLY SET TO 

A NON- ZERO VALUE 

ANGLE: BA - BASIC LOAD VAL=0.0 OR 1.0 

A - ANGLE OF ATTACK (a) IN DEGREES 

B - ANGLE OF SIDESLIP (g) IN DEGREES 

P - ROLL RATE IN RAD/SEC 

Q - PITCH RATE IN RAD/SEC 

R - YAW RATE IN RAD/SEC 

T - THICKNESS VAL=0.0 OR 1.0 

DA# - ANTISYMMETRIC CONTROL SURFACE DEFLECTION NUMBER # 
DS# SYMMETRIC CONTROL SURFACE DEFLECTION NUMBER # 

S: SEC1 TO SEC2 SECTIONS OR COMPONENTS 

M: MACH NUMBER 

SA: ABSCISSA SCALING FACTOR 

0.0 FREE SCALING. DEFAULT 0.10 

SO: ORDINATE SCALING FACTOR 

0.0 FREE SCALING. DEFAULT 0.0 
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ABBREVIATIONS: 


NONE 


FUNCTION: 


SET UP FOR PLOTS OF QUANTITIES WHICH CAN BE ILLUSTRATED 
AS CHORDWISE AND SPANWISE DISTRIBUTIONS. 

SETS OF ANGLES AND SECTIONS CAN BE PRE-SET USING ADD AND 
S .THEN PLOTS CAN BE SELECTED FOR DIFFERENT COM’S AND 
MACH NUMBER. 
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-PURGE- 


FUNCTION PURGE 



OPERANDS : NONE 


ABBREVIATION: PURGE 


FUNCTION: RELEASES ALL COMPONENTS IN THE LOCAL FILE. 
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- RENAME - 


. FUNCTION RENAME 


RENAME, COMP[REC] [,FILE] 


OPERANDS: COMP - COMPONENT NUMBER OR RECORD NUMBER REC 

FILE - PERM [p] , {DEFAULT FILE J LOCAL £ L .] 

ABBREVIATION: RENA 

FUNCTION: CHANGE COMPONENT NAME 

ADDITIONAL RESPONSES :OLD NAME: "COMPONENT NAME" ENTER NEW NAME- 

NEW NAME. 
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-RENUMBER- 


FUNCTION RENUMBER 



OPERANDS: 


OLD CQMPN: OLD COMPONENT NUMBER, OR RECORD NUMBER [REG] 
NEW CCMPN: NEWCOMPQNENT NUMBER . 

FILE: j^PERM [p] , [DEFAULT FILE]J JlOCAL [l]J 


ABBREVIATION: RENU 


FUNCTION: 


CHANGE THE COMPONENT NUMBER. 
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-SAVE- 


FUNCTION SAVE 


OPERANDS: 


ABBREVIATION: 


FUNCTION: 


save[,file] 


FILE - PERM [P] [DEFAULT FILE] 
LOCAL [L] 


SAVE 


RETURN COMPONENT IN CORE TO ITS RECORD LOCATION IN 
REQUESTED FILE. IF NO COMPONENT WITH THE SAME NUMBER 
IS FOUND. THE COMPONENT WILL BE CATALOGED. 
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^ SLENDER - 


FUNCTION SLENDER 


SLENDER, COO. [RECl] .... ,CCMN(RECN] [FILE] 


OPERANDS: CCM1: 

FILE: 

ABBREVIATION: SLEN 

FUNCTION: CONVERTS UP TO 6 SELECTED COMPONENTS INTO A SINGLE 

SLENDER BODY FOR ISOLATED BODY CALCULATIONS. THIS 
FUNCTION MUST BE USED TO CREATE COMPONENTS FOR THE 
ISOLATED BODY SOLUTION. 

ADDITIONAL RESPONSES 

* SLENDER BODY SIMULATION 

* ENTER COMPONENT NAME -"SLENDER BODY" 

* INPUT COMPONENT NUMBER -"700.0" 

* COMPONENT PARAMETERS FOR SLENDER BODY COMBINATION 

* ENTER FACTOR,NPASS -*0.7,2" 


COMPONENTS TO BE BLENDED INTO SLENDER BODY GEOMETRY. 
ONE OR MORE MUST APPEAR 
[RECl] RECORD NUMBER MAYBE USED INSTEAD OF 
COMPONENT NUMBER 

FpERM [P] , DEFAULT! , f LOCAL [L]l 
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-SLENDER- 


FACTOR: FRACTION OF DISTANCE TO NEXT POINT ON PRESENT CROSS 

SECTION TO BE COMPARED AGAINST DISTANCE TO POINTS ON 
OTHER CROSS SECTIONS WHEN GENERATING OUTLINES. ONLY 
NECESSARY WHEN TWO COMBINING COMPONENTS ARE CUT AT 
THE SAME CROSS SECTION 


NPASS: NUMBER OF POINTS ON CROSS SECTION THAT SHOULD BE 

SKIPPED AFTER OUTLINE PASSES FROM ONE COMPONENT 
TO ANOTHER. 
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-STATUS- 


FUNCTION STATUS 


OPERANDS: 


ABBREVIATION: 


FUNCTION: 


STATUS 


NONE 


STAT 


RETURNS THE NUMBER AND NAME OF THE COMPONENT IN CORE 


102 




-TERMINAL- 


FUNCTION TERMINAL 


TERMINAL 


OPERANDS : NONE 


ABBREVIATION: TERM 


FUNCTION: USING KEYS GIVEN BY THE PROGRAM, THE OPERATOR INPUTS 

SIMPLE COMPONENT GEOMETRIES CONSISTING OF 
GEOMETRICALLY SIMILAR SECTIONS. 
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-TERMINAL- 


ADDITIONAL RESPONSES AND INPUT: 


* 1 FPLAN * 2 HPLAN * 3 FELIP * +- 4 HELIP * 5 RECT * + -6 TRIANG * 7 END - 


THE USER INPUTS THE NUMBER INDICATING THE TYPE OF COMPONENT FORM DESIRED: 

1 : TYPE 4 FULL SURFACE 

2 : TYPE 3 HALF SURFACE 

3 : TYPE .1 OR 2. FULL ELLIPSE 

4 : TYPE 1 OR 2 HALF ELLIPSE + 0-0 

5 : TYPE 1 OR 2 RECTANGLE 

6 : TYPE 1 OR 2 TRIANGLE +A -V 

7 : END 


* INPUT COMPONENT NAME: "FUSELAGE” (16 CHARACTERS MAXIMUM) 

* INPUT COMPONENT NUMBER: "100* (>50) 

SURFACE COMPONENT INPUT : (EXAMPLE SHOWN IS FOR A TYPE 4 INPUT) 

* FULL PLANAR: X(CBAR/4) ,YO(ROCT) ,Z0 (ROOT) -”250, ,30,0." 

* AR,S, TAPER, SWEEP (DEG) ,DIH (DEG) -"4. 5, 250. ,0.32,45, 45.0" 

* IAF (1,2, 3, 4, OR 5), (T/C)IN, (T/C)0UT-"2, 0.04, 0.05" 
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WHERE IAF VALUES ARE: 


1 : 6SA0XX 

2 : 64A0XX 

3 : SUPERCRITICAL 

4 : HEX - AIRFOIL 

5 : Bl- CONVEX AIRFOIL 

BODY COMPONENT INPUT: (EXAMPLE IS A CENTERLINE ELLIPSE) 

* NON PLANAR COMPONENT XRO, YRO, ZRO -"0.,0,0.“ 

* NCS.X &) ,A(1) ,A/B(1) , . , . ,X(NCS) ,A(NCS) ,A/B(NCS) 

-'4,0. ,0. ,1. ,200. ,1000. ,0.6,400. ,2500. ,0.8,600. ,1500. ,1.5* 

MAXIMUM HEIGHT ABOVE COMP. REF. LINE: 36.4183 

MAXIMUM HALF BREATH OF COMPONENT : 21.8510 
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-TIME- 


FUNCTION TIME 


ABBREVIATION: 


FUNCTION: 


TIME 


NONE 


GIVES WALL TIME, CPU TIME, AND APPROX. COST SINCE START 
AND SINCE LAST CALL. 
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-TITLE- 


FUNCTION TITLE 


TITLE [,FILE] 


OPERANDS: FILE - PERM [P] [DEFAULT FILE] 

LOCAL [L] 


ABBREVIATION: TITL 


FUNCTION: INPUT OR CHANGE TITLE OF REQUESTED FILE 


ADDITIONAL RESPONSES: OLD TITLE: TITLE 

ENTER NEW TITLE 

TITLE ,[P=L] ,[L=P] 


P - L - PERM TITLE = LOCAL TITLE 
L = P - LOCAL TITLE = PERM TITLE 
(space) (CR) - RETAIN PRESENT TITLE 
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-UNITS- 


FUNCTION UNITS 



OPERANDS: FILE: PERM[p] [DEFAULT FILE] 

LOCAL [L] 

ABBREVIATION: UNIT 

ADDITIONAL INPUTS: 

UNITS 1 -METERS 2 -ENGLISH 3-CHOTMETERS - "UNI" 
UNI - NUMBER INDICATING UNITS DESIRED 
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GLOSSARY 


CARD: Data description using card input. Geometric definitions are 

similar to those used in TERMINAL. 

Component: A geometric description of a vehicle or portion of a vehicle 

which by itself can be considered a complete unit. 

Cross-Section: A cut through a component made perpendicular to the 

principle axis of the component . Planar components have a principle axis 
defined by the leading edge trace in the Y-Z plane. Non-planar components 
have a principle axis defined in the x direction. 

Digitize: To define components with cross-sections of the component using 

a graphics tab let. and an appropriate input device. 

Interference: Refers to a paneled component used in conjunction with 

isolated body, i.e. interference shell. Interference shells default to 
surface components when an isolated body is not present. 

Body: Refers to components such as a fuselage, nacelles, etc. 

Surface: Refers to components such as wing, tail, etc. 


Slender: Refers to a component constructed from a single or set of body 
components using the function SLENDER. 

Terminal: Refers to components constructed by inputting simple geometry 
definitions such as cross-sectional and width/height for body components 
and for surface components. 
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APPENDIX I 


DATA DESCRIPTION 
FOR CARD INPUT GEOMETRY 
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NUMBER 



IDENTIFICATION DESCRIPTION DO NOT KEY PUNCH 


TITLE (72 Characters) 



NC - No. of Components (< 50) 


COMPONENT NUMBER 1 



COMPONENT NUMBER 50 


(Last card before first component data must have a minus sign 


in column 1 


(Last card of each component must have a minus sign in column 1) 




IDENTIFICATION DESCRIPTION 


NOT KEY PUNCH 



COMPONENT NUMBER (> 50.0) 


COMPONENT NAME (16 Characters) 


BODY COMPONENTS 



NCS Number of cross section ( ^18 ) 


Relative origin of component 


1:C.L. Full ellipse, 2:off-set full ellipse 


RLC 

(Locations 105-199 are for section data) 


Al 

cross sectional area 


Dl 

h/2 (full height of A and O ) 




o . o 


0 . 0 


Section origin relative to X0,Y0,Z0 



7 M I . 0 


0 . 0 


(o.2. 


I 90 



RLC + NCS* 5 


ANCS 


dncs 


*R 


Yr 


ZR 

























DESCRIPTION DO NOT KEY PUNCH 


COMPONENT NUMBER (> 50.0) 
COMPONENT NAME (16 Characters) 


SURFACE COMPONENTS 


NCS Number of cross sections ^ 17 sections 

Xg 

Yq Relative origin of component 

Zcl 

Ctype 4.0 wing [full surface]; 3.0 vert, [half surface] 


k 1.0-65A, 2.0-64A, 3. 0- supercritical t 

AIRFOIL TYPE 4.0-HEX, 5.0 Bl-CONVEX, 0,.0-USER INPUT Z/C'S 

UPPER SURFACE Flag; if (DA(105) = 0.0 set = 1.0 

LOWER SURFACE Flag; if syirrnetric leave 0.0 1 ■' 

NO. X/C'S INPUT; if DA(105) =0.0. 


AR 

S 

TAPER 

SWEEP 


DIHEDRAL 
















IDENTIFICATION DESCRIPTION DO NOT KEY PUNCH 



RLC 115 - 199 are for section data 


T/Q4AX , x 


CHORD 


CONSECUTIVELY, ONE 


XR l card PER CROSS SECTION 


Yr Coord, of L.E. relative to Xo,YoZ 0 ]2ND SECTION IS CARD 


ZR * [l20, 3RD 125, ETC. 


RLC 200 - 219 are for x/c‘s 


X/C X 





DA (220- 399) Are for upper and lower surface Z/C’s. The 


Z/C's of one section are followed directly by the Z/C's for the 


lower surface. 


(fill in only if UPPER SURFACE flag is >0.0) 
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APPENDIX II 


PROGRAM COMPILATION, LOADING AND SIZE 

The Preliminary Aerodynamic Analysis System is comprised of three major 
programs. The input/ output (GEOMETRY) routine which is interactive and 
the slender body (SLENDER) and lifting surface (LIFTING) programs which are 
designed to run batch. The names in parenthesis refer to the CDC file names 
which will be used to identify the source codes of the three programs in what 
follows,. 

GEOMETRY is loaded with a segmented load which reduces the core 
requirements for the loaded program to approximately 7 OK octal on the CDC 175 
system. In addition, the coding itself must be linked with libraries 
containing TEKTRONIX PLOT- 10 and GRAPHICS TABLET enhanced graphics software. 
Listing II -1 illustrates the BATCH submitted job control listing (JCL) to 
compile and link the GEOMETRY source code into an executable load deck 
which can be run on a TEKTRONIX 4014 graphics teiminal. Listing II -2 
illustrates the segmentation listing (SEGMENT) used in the link. The 
resulting load deck (listing II -1) is cataloged as LANGLOAD and is now 
available for interactive execution. 

LIFTING is loaded with a segmented load to reduce the core needed to 
run the program in batch. The JCL needed to compile and link LIFTING is 
shown in listing II-3. The segmentation list is shown in listing II-4. The 
resulting load deck is filed under LIFTLOAD for later batch execution. 
LIFTLOAD has an operating core size of approximately 114 K octal. 

SLENDER is loaded directly by compiling and linking the program, see 
listing II -5. The resulting load deck is then stored as SLENLOAD for later 
batch execution. SLENLOAD has a core size of approximately 100 K octal. 
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100=PEDAJ,T120, 10200, CM120000,P4. 

110=ACCT (DIVAN 0018422401*0111013) 
120=REQU£ST,PAD1 , *PF. 

130=ATTACH, ONE, LANGLEY, ID=D0500. 

14 0=FTN , I=ONE , B-LOAl ,L=0. 

150=MAP,PART. 

16 0=ATTACH , PLOT10 , PLOT10LIB , ID=APLOT10 , SN=SYSAPP . 
170=LIBRARY, PLOT10 . 

180=ATTACH , SEG , SEGMENT , ID-D0500 . 
190=SEGLOAD(I=SEG,B*PAD1,LO=DT) 

200«LOAD(LOA1) . 

210=sNOGO* 

2 2 0 =C ATALOG , PAD1 , LANGLOAD , ID=D0500 , RP*10 0 . 
240=DISPOSE , OUTPUT , PR*I AJ . 

250=EXIT. 

26 0 as DISPOSE , OUTPUT ,PR-IAJ. 


Listing II -I JOB CONTROL LISTING for Compilation 
and Linking of Geometry program, 
GEOMETRY 
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100=START TREE MAIM 
110=MAIN INCLUDE MAIN 
120= LEVEL 

130=EXEC TREE GEOEX 
140=GEOEX INCLUDE GEOEX 
150=DIGITI TREE DIGIT 

1 6 0=DIGIT INCLUDE DIGIT ,ORIG , TABOFF , TABARM , MULPNT , TABOFF , GETPNT 
170=SAVE TREE SAVEB 
180=SAVEB INCLUDE SAVEB 

190=VIEW TREE VIEW3- (INTERF, TRANS, INSTR2, AREA) 

200=VIEW3 INCLUDE VIEW3 
201=INTERF INCLUDE INTERF 
20 2= TRANS INCLUDE TRANS 

203=INSTR2 INCLUDE INSTR2 
204=AREA INCLUDE AREA, LIMIT, SHIFT 
210=VIEW3 GLOBAL V3 

220=EDIT TREE SECT- (EDITP,EDITN, LIST) 

230=SECT INCLUDE SECT 
240=SECT GLOBAL ED 
250=EDITP INCLUDE EDITP,PLOTOT 
260=EDITN INCLUDE EDITN 
270=LIST INCLUDE LIST,EDPUT 
280=DISPL TREE DISP 

290=DISP INCLUDE DISP, DISPS , INFO, GREEK 
300=SLENDR TREE SLEND 
310=SLEND INCLUDE SLEND 
320=TERM TREE ICON 

3 3 0=ICON INCLUDE ICON , EMLORD , MATRIX , TRAP , XYZ , WCAD , DECRD 
340=ANALY TREE ANAL 
350=ANAL INCLUDE ANAL 

36 0=ANAL GLOBAL PLOT , CDAT , CFIX ,CFLO , CLIM-S AVE 
361 =GR APH TRE E PLOTT- ( LOADS , FORC ) 

362=PLOTT GLOBAL PLOTT-SAVE 
380=LOADS TREE CPPR- (PLEXC ,CPWING) 

390=CPPR INCLUDE CPPR 

400=PLEXC INCLUDE PLEXC 

410=CPWING INCLUDE CPWING, LABEL 

420=FORC TREE FORCE- (TOTAL, SHEET, FOREX) 

425=FORCE INCLUDE FORCE 
430=TOTAL INCLUDE TOTAL, DRGPOL 
440= LEVEL 

450=LOOK TREE LOOKUP- (SPLINE, FITT) 

460=LOOKUP INCLUDE LOOKUP, FIT 
470=SPLINE INCLUDE SPLINE 
480=FITT TREE HSHLDR-FITWT 
490=HSHLDR INCLUDE HSHLDR 
500=FITWT INCLUDE FITWT 
510=DRAG TREE CPLOT- ( WAVE , VISCO) 

520=CPLOT INCLUDE CPLOT 
530=CPLOT GLOBAL CXXX 

Listing II -2 Segmentation Listing of Permanent File 
SEGMENT (See Listing II-l line 180) 

Used to Link GEOMETRY Program 
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540=WAVE 

5 50=WDRAG 

560=WDRAG 

570=DQWS 

580=XYZR 

590-AREAW 

600=XACALW 

610=VISCO 

620=VDRG 

6 30= VDRG 

640=AREAS 

650*WSVP 

66 0*1 DATA 

670=TLU 

680=LIFT 

690=LSET 

700=PLOT 

710=PLOTIT 

712=EMLOR 

714=EMLO 

720 = 


TREE WDRAG- (DQWS,XYZR- ( AREAW , XACALW) ) 

INCLUDE WDRAG 

GLOBAL CTHE ,COPT 

INCLUDE DQWS,SDQ,COD 

INCLUDE XYZR, LIMITW 

INCLUDE AREAW, SHIFTW 

INCLUDE XACALW 

TREE VDRG- ( AREAS , WSVP , I DATA ,TLU) 

INCLUDE VDRG 

GLOBAL COUT,CLUN 

INCLUDE AREAS , LIMITS , SHIFTS , WET 

INCLUDE WSVP 

INCLUDE IDATA 

INCLUDE TLU 

TREE LSET 

INCLUDE LSET, SMOOTH 

TREE PLOTIT 

INCLUDE PLOTIT, GRID, SCALIT,MAXIT, CODIM 
TREE EMLO 
INCLUDE EMLO,ISIMEQ 
END 


Listing II -2 Completed. Segmentation Listing of 
Permanent File SEGMENT (See Listing 
II-l line 180) Used to Link GEOMETRY 
Program 
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100=PEDAJ,T120, 10200 ,CM240300,P4. 

110*ACCT (DIVAN 0018422401*0111013) 

1 20=REQUEST , UDPL , *PF 
1 30=ATTACH , ONE , UDPFORT , ID=D0 5 0 0 . 

1 4 0 =FTN , I ONE , B =L0A 1 , L= 0 . 

150-MAP, PART. 

16 ATTACH , SEG,SEGUDP, ID«D0500 . 
170-SEGLOAD (I*SEG,B=UDP1,L0=DT) 

180 =L0AD ( L0 Al ) . 

190=NOGO. 

20 0-CATALOG, UDP1 , UDPLOAD, ID=D0 5 00, RP=100. 
210=DISPGSE, OUTPUT, PR=IAJ. 

220*EXIT. 

230 =DISPOS E , OUTPUT , PR= I A J . 


Listing II -3 JOB CONTROL LISTING for Compilation 
and Linking of Lifting Surface 
Program, UDPFORT 


123 



100=START 

TREE 

UDP 

110=UDP 

INCLUDE 

UDP , PAGE , CODIM , FI LL 

120= 

LEVEL 


130=MAT 

TREE 

MATRIX- (VORTEX, SOURCE) 

140=MATRIX 

INCLUDE 

MATRIX, BMTX 

150=MATRIX 

GLOBAL 

Will -SAVE 

160=VORTEX 

INCLUDE 

VORTEX 

170=SOURCE 

INCLUDE 

SOURCE 

180=BOD 

TREE 

BOUND 

190=BOUND 

INCLUDE 

BOUND 

200=JT 

TREE 

JET 

210=JET 

INCLUDE 

JET 

220=SOL 

TREE 

SOLVE 

230=SOLVE 

INCLUDE 

SOLVE , MSOL , REVERS 

240=FOR 

TREE 

FORCE- ( PRESS , INTEG , DRAG- (TRIM , VTX ) ) 

250=FORCE 

INCLUDE 

FORCE 

260=PRESS 

INCLUDE 

PRESS 

270=INTEG 

INCLUDE 

INTEG, BDYLD,JETLD 

280=DRAG 

INCLUDE 

DRAG 

290=TRIM 

INCLUDE 

TRIM 

300=VTX 

INCLUDE 

VTX, NORMAL 

310=WAV 

TREE 

WAVE 

320=WAVE 

INCLUDE 

WAVE 

330=PRI 

TREE 

PRINT 

340=PRINT 

INCLUDE 

PRINT 


Listing II -4 Listing of Permanent File 
SEGUDP (See Listing II- 3, 
line 160) Used to link 
UDPFORT Program 
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100=PEDAK / T120, 10200 f CM120000 , P4 . 

110=ACCT (DIVAN 0018422401*0111013) 

1 20=REQUEST, SLEN / *PF . 

130=ATTACH /ONE/ SLENDER/ ID=D05 00 • 

140=FTN/ I=ONE ,B=L0A1 /R*0 . 

150 =M AP / PART • 

160=LOAD (LOAl ) . 

170=NOGO (SLEN) . 

1 8 0=CATALOG , SLEN , S LENLOAD , I D=D0 5 0 0 r RP= 999. 
190=DISPOSE /OUTPUT/ PR=IAJ. 

200=EXIT. 

210=DISPOSE /OUTPUT /PR=IAJ. 


Listing II -5 JOB CONTROL LISTING for 

Compilation and Linking of 
Slender Body Program, 
SLENDER 
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APPENDIX III 


SAMPLE SESSIONS OF SYSTEM OPERATION 

The four sessions presented illustrate typical program operations 
involved in aircraft configuration input, editing, analysis and display of 
results. There are also examples on how to store and retrieve files created 
using the programs of this system. 

The figure on the following page illustrates the type of drawing 
required for a configuration. A side view and plan view are necessary along 
with cross section cuts at various axial station, drawn to the same scale as 
the top and side views. These sections are digitized to simulate the vehicle 
fuselage. 

Of the four sessions which follow, the first three are for the fighter 
of Figure III -I. The last session is for a small trainer design to 
illustrate typical lifting surface displays. 
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SESSION I 


CONFIGURATION INPUT AND EDITING 


CONNAND- ATTACH, LANO, LANOLOAD, ID-DtfM . 

PF CVCLC NO* • Ml 
CONN AND- CONNECT, INPUT, OUTPUT. 

CONNAND- ■ ETL, 18# 

CONNAND- SCREEN, 13* 

COflNAND- REWIND, LANQ 

CONNAND- LANO 

* INPUT TERN INAL CHARACTER SPEED l 3* 


To initiate the input /output program (LANGLOAD) the 
user first attaches the machine language file (COMMAND 1). 

The input and output files must be connected to the teminal 
(OTJviAND 2). If the system has a default execution time 
limit, the user should over ride that (COMMAND 3). To take 
advantage of full screen writing the user may have to set the 
screen size (COMMAND 4) . The user then rewinds the machine 
file (CONMAND 5) and executes the program (COMMAND 6) . 

The program will then request the transmission speed of the 
terminal in characters per second (bottom line). 
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t KEY * X0,XRAX t 

• 400 

PEh * PLANFORR XO, PROFILE XO, PROFILE XRAX * 
AXIS OK t 1 YES * 3 HO 

*» INPUT |AV PARARETER0 

? KEV t TYPEU,8,3 OR 4> # NSE0,NOUT<n,..,NOUT<N$EO> 

INPUT CORP NURIERl 

100.1 

INPUT CORPONENT NARE* 

ifi futt Ugt 


D 


D 


In the digitize function the program requests a scaleable 
length, indicated by the input of XO,XMAX as the limits 
on the fuselage reference axis. These limits need not 
encompass all cross sections to be digitized but can be 
any two convenient locations on the reference line . 

The program now asks the user to digitize the (XD, y * 0) 
point of the plan view, (TOP), the (XO,Z = 0) point of the 
profile view (SIDE) and the (XMAX, Z * 0) point of the 
profile view. The routine will display the three input 
points and the user can check them for consistency and 
orthogonality. 

The program now requests input of the component type, the 
number of segments per cross section (NSEG) , and the wetted 
flag for each segment (NOUT = 1 wetted, -1 un- wetted) . 

In the example a type 1 (centerline body) with 5 segments; 
segments 1, 2 and 4 wetted, 3 and 5 unwetted has been 
selected. 

The program now requests input of the component number and 
name. 


P 
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The user digitizes points around each cross section, 
entering S(CR) after each completed cross section.** 

In this case the first try at the cross section resulted 
in an incorrect number of input segments (denoted by double 
points digitized by user) and the diognostic seen between the 
two horizontal arrows was printed and the tablet was turned 
on automatically so the section could be re- digitized. 

The input was corect this time and a cross section display 
(left side of screen) and tabulation of input data 
(upper right of screen) are displayed on the tube. 


The user now responds to the lower left hand instruction; 
in this case and in the next figure the option to digitize. 
The next section is selected. In the third figure of this 
set the user selects the end option or (3) 

** The digitizing routine the user may supply may not 
contain the S(CR) to end digitizing and the user’s act- 
ions may need to be adjusted accordingly. 
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t SMOOTH OPT I ON I 0-NO, ENTER ORDER-VES * 
7 

0 

ttOKft 

catalog 


COMPONENT 100.00 ALREADY EXISTS 

ENTER NEU NUNIER OR 0.0 TO RE-USE OLD LOCATION- 

? 

700. 

noKtt 

tiaa 


UALL TINE I 
DELTA UALLi 
«OKM 
f lias 

FILE* PERN 
REC COMP NO. 


14.70 CPU TINE i 
0.82 DELTA CPU! 


NAME 


1 

8 

3 

4 

5 

6 
7 
I 
9 

10 

ttOKtt 

rflfi 


109.00 FOR FUSELAGE 

101.00 AFT FUSE 

103.00 NACELLE 

800.00 IN9D UING 

801.00 OUTID UXNG 

300.00 LOUER TIP FIN 

301.00 UPPER TIP FIN 

400.00 LOUER VERT 

401.00 UPPER UERT 

500.00 VERT PYLON 


0.98 

0.01 


TYPE 


3 

8 


The program inquires whether the user wants to smooth 
the input points of each segment or simply spline fit. 
In this instance the user has selected the spline fit 
only (0) 

The program returns to the **OK** mode and the user 
indicated he wishes to catalog the new component. The 
program responds that the component already exists and 
the user enters a new number (700,0) to avoid writing 
over the old location. 

The user requests the time for this session, then asks 
for a listing of the files in the permanent file. The 
user then requests the digitizer again. 
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KEV # TOC (SET NEGATIVE IF REFERENCE) 

PEN X IE (PROFILE), IE AND TE (PLANFORN CHORD LENGTH) 
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The program asks for the airfoil thickness. If TOC 
is less than zero the program expects a digitized 
airfoil. Since the user has input a positive thickness 
the program assumes the previous input or, in this 
instance, the default (64A0XX) airfoil as the desired 
section . 

The program turns on the tablet and the user inputs the 
Z location of the section (point 1) , the Y and X of the 
leading edge (point 2) and the Y of the trailing edge 
(point 3). 


The scaled section is then drawn in to show that the 
points have been input and the section characteristics are 
shown in the upper right hand comer. 



\ 

D 

> 1 NEXT SECTION t 8 REPEAT SECTION 8 3 END COUP. 


Note : if a negative thickness is selected the user would 
then input a cross-section starting with the reference 
chord length then the cross-section itself. Since the 
program calculates local dihedral and section twist about 
the first input point on the wing, the first point should 
lie on the reference chord line. 
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SMOOTH OPTION* l-NO, ENTER ORDER-VES t 

L.E. RADIUS l METHOD Hit METHOD 8t 8 t OFF f 3 t 


The user now selects to smooth the input sections, in 
this case to better calculate the leading edge, by inputting 
a "3” for third order smoothing. 

Airfoil section can be fit with a leading edge radius 
polynomial built into the smoothing routine. In this 
instance the user selected option 2 which allows the 
routine to determine the leading edge radius. 

The next two figures present the display showing the 
smooth fit and the spline fit which the user can select 
from. (1, or 0). In this case the user selected the smooth 
in all instances and then went directly into edit (lower 
left and top middle of second figure) on the new component . 
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n C0NF10UNATI0N CONCEPT A/A FIOHTEA CANAW 

* OPTIONS! NO CHANGE • I CHOASUISE FIT SECTIONS I t CHANGE SPAN STATIONS * * 
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tt CONFIGURATION CONCEPT A/A FIGHTER CANARD 

t OPTIONS* NO CHANGE i * CHORDUISE FIT SECTIONS 1 t CHANGE SPAN STATIONS 2 f 



INSTRUCTIONS 

POSITION SPAN STATION 
I* INTERPOLATE STATION 
Dl DROP STATION 
E* CHORD ELONGATION 
LOCATE L.l. 

LOCATE T.E. 

SPACE* END STATION INPUT 
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ft 


CONFIGURATION CONCEPT A/A FIGHTER 


ft OK ft 

catalog 

ttOKtt 

rfUj)*tkraa»all 


In this final figure of the series the program 
has displayed a view showing the final work 
performed on the component and returns the 
users to the **OK** mode. The user catalogues 
the component then asks for a three view 
display of all the components . 




CONFIGURATION CONCEPT A/A FIGHTER 



The user elects to display cross-section 
cuts of the configuration by selecting option 
two. By using the vertical graphics cursor to ^ 
locate the x- station and the horizontal to position 
the vertical screen location of the displayed section 
cut, and working from left to right to avoid number 
overlap, the results shown in the next figure are 
obtained. 
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rtOKtt 
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FItCi 


PERN 


REC 

COUP NO. 

NAflE 

1 

1H.N 

FOR FUSELAGE 

2 

lei.ss 

AFT FUSE 

3 

183.88 

NACELLE 

4 

888 .88 

INID UINQ 

5 

881.88 

OUT ID UINQ 

6 

388.88 

LOUER TIP FIN 

7 

381.88 

UPPER TIP FIN 

8 

488.88 

LOUER UERT 

9 

481.88 

UPPER UERT 

tt 

$88.88 

UERT PYLON 

11 

ItOKtt 

688.88 

CANARD 

lnitrftranc* l 8 888*688 


TYPE 

1 
i 
1 
4 
4 
3 
3 
3 
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The program returns to **OK** and the user asks for a list of 
components (files). 

He now requests four components to be used in making an interference 
shell, records 7 and 2 and components 200 and 600 (Note the mixing 
of records and components) 

In the next figure the user has been supplied with a three view with 
bodies dashed in as references only. Surfaces in this view can be used 
to line up or "connect” construction lines of the inteference shell to. 

The user first locates X leading edge and Y in the planview and the 
number of the connecting or alignment component (numbers displayed 
at centoid of surface components, or a "0" for a free line. 

The user next locates X trailing edge and Z and presses the type of 
surface he is constructing, (I, E, or J) the program will than display 
the constructed span line in the two views . 

On subsequant passes to construct additonal lines, the user need only 
locate Y and Z if an interference shell is selected since X leading 
and trailing edge are fixed by the first line. With the standard 
boundary condition the X's must always be input and with a jet flap 
the leading edge X must match a connecting body the trailing edge X 
is a function of a user input length of the connecting chord. 

The resulting line construction is shown in the following figure, a 
space used in the first step ending the routine (see menu in middle 
right side of page). 



iir »r 

*qmoo 'SHnoi awMMUf »i 

33M3)BJNliMI 1 1 

isau tc 

3M3 ONniVHl 3XV301 tl 
NOIXVOOl 2 3XV001 tl 
(AHA MID It 2li 

sxnduoo mnmm »iow« 

3NI1 3Ui it 
XM3MMM03 0MIX33NM03 It 

SS3)M K 
joas Miatm 3xvooi n 
N0IXV301 A 1XV301 It 
(faXAIWID I 43X1 


SMoixontaSMi 



•*MXUtifUnN XHBNOdUOO s 
n«M« nw? 40)1 3UVN XM3WOdUOO t 





M3XH9IJ W/V Xd33KI00 NOUMflOIiNOO 



CONFIGURATION CONCEPT A/A FIGHTER 




INSTRUCTIONS 

. .Jiff 1 (PLANVIEW) 

It LOCATE Y LOCATION 
1* LOCATE LEADING EDGE 
ji PRESS 

ft CONNECTING CONPONENT 
•i FREE LINE 

SPACE l PANELING CONPLETE 

STEP II (SIDE VIEW) 

It LOCATE S LOCATION 
It LOCATE TRAILING EDQE 
3» PRESS 

It INTERFERENCE PANEL 
Et STANDARD SOUND. COND. 
Jt JET FLAP 


ft CONFIGURATION CONCEPT A/A FIGHTER 


FOR INTE SHEtt ft 


in 

O 



I NO CHANGE • « CHOROUISE CONSTRAINTS 
* FULL COS SPACING 4 f 
t ENTER NO. OF CHOUSE PANELS * 
t PRESENT NUNIER* 1 t 
? 

S 


1 S EVEN SPACING 2 t HA 


The user goes into edit for the interference shell 
selecting even spacing for the shell which 
comes with one panel from construction. 


The display of the final view is shown in the 
following figure with the five evenly spaced panels 
inserted. 


The user then catalogues the shell and asks to 
construct the aft shell, which is similar and 
not shown here. 

The final view for this shell is shown in the 
following figure. The command input at the middle 
of this figure is to display the interference shells 
and the wing and canard to illustrate the connections 
between surface components. This is seen in the last 
and final figure in this series. 
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tf CONFIGURATION CONCERT A/A FIGHTER 


*AFT INTE SHELL ** 


ttOKft 

ciUUf 

tfOKtf 

nut 

FILE* PERN 
REC COUP NO. 

NANE 

TYPE 

1 

199.99 

FOR FUSELAGE 

1 

2 

191.99 

AFT FUSE 

1 

3 

193.99 

NACELLE 

1 

4 

299.99 

INBD WING 

4 

6 

291.99 

OUTID WING 

4 

6 • 

399.99 

LOUER TIP FIN 

3 

7 

391.99 

UPPER TIP FIN 

3 

8 

499.99 

LOUER UERT 

3 

9 

491.99 

UPPER UERT 

3 

1# 

S99.99 

UERT PYLON 

2 

It 

699.99 

CANARD 

4 

12 

799.99 

FOR INTE SHELL 

S 

13 

792.99 

AFT INTE SHELL 

s 

ttOKXt 

dlsp,t,299,291. 

11,12*13 





CONFIGURATION CONCEPT A/A FIGHTER 



1W.7 345*5 


I TRANSLATE 1 t CROSS SECTIONS 8 * END • t 
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NORNAL JOB COMPLETION 

CONN AND- catalog,tapal,gaoaatry,ld*ttit»i,rp»999,c*«opan,ax«rila. 
NEU CYCLE CATALOG, CY»Dtl 

CONN AND- cats log, tapa3,plotf i la, ld*tMtM,rp*999,cn»opan,ax*f 1 la. 

NEU CYCLE CATALOG, CV-M1 
CONNAMD* logout 


In this frame the user has exited from the system 
and is cataloging his geometry and plot files for 
later use. The use of mass storage files 
necessitates storing the files so they can be 
extended if necessary on future runs. The addition 
of 

Qsj = open, ex=file 

to the end of the catalog statement assures that this 
will occur. 



SESSIQU II 

WAVE DRAG AND VISCOUS DRAG ANALYSIS 


CONflftNO- itiach« tang. langtead. 

PF CYCLE NO. -Ml 
COMMAND- connect. input. output • 

COMMAND- attach, tape 1 . geeeetr y •itf*titttt,pw • e pan .file. 

PF CYCLE NO. -Ml 

COMMAND- attach. taped, p lot f l la. id-tttMt, pw-open. file. 

PF CVCU NO. -Ml 

COMMAND- attach. taped, wevedr eg. id I OPTIONAL ENTRY SEE FOOTNOTE! 

PF CYCLE NO. -Ml 

The user attaches the program load file, the geometry file 
and plot file from the previous session. If a lifting surface 
solution had been performed the user might also want to attach 
the file containing pressures for the wave drag due to lift 
solution, however this is optional. 


COMMAND- rewind.leng. 

COf*IANO> tang. 

* INPUT TERMINAL CHARACTER SPEED I 3t 


(FILE UAUEDRAG (TAPE* ) 
UAWE DRAG DUE TO LIFT 


ONLV NEEDED UHEN A UDP SOLUTION HAS BEEN OBTAINED AND 
IS REQUIRED FRON ANALVSIS) 
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DISPOSITION OF PEW file (NEU/OLD)! 
old 

DISPOSITION OF PLOT FILE lN£U/OLD)i 
old 

MOKtt 

HIM 

OLD TITLE* 



The user informs the program that both the perm file 
(geometry) and plot file already exist and to read in their 
contents record. 

The title of the perm file is requested and the local file 
title is set equal to the perm file title (£ = p). 

Tiles in the perm file are listed and components in records 
1-11 are attached. Note: an easier way would be to attach 

all components and delete the last two from the local file. 
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The keyboard ANALY** indicates the analysis function and the 
user inputs the Mach numbers to be used in this exercise. 
Note the use of Mach 0.60 for wave drag to include the 
empirical drag divergence interpolation in the transonic 
regime . 

The user requests wave drag and viscous drag solutions, 
executes these solutions by entering "calculate" or "C". 

The program displays the reference data which it has stored 
giving the user an opportunity to change the data, then 
enters the first analysis routine, viscous drag. 

The scale requested is based on vehicle inches/drawing inch, 
of digitized configuration. A modem fighter for example 
would be 20 scale (20 inches/ inch) , the value entered here. 

The use of the word plot here refers to seeing selected 
cross section cuts to remove perimeter segments that are 
unwetted. In most cases the third option (51 cuts) gives 
an accurate estimation of wetted area for viscous drag 
estimates. 
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The display shown in the next two figures are the perimeter versus X- station and the cross sectional 
area versus X- station. The scales listed for each plot are valid except for the top horizontal line and 
right vertical line and were selected based on a fighter sized vehicle. 

The third figure in this series shows the component data necessary to run viscous drag. The skin roughness 
is estimated with a sand grain height (in feet). The pressure, temperature and, Mach (disp) (the Mach 
number for which a component drag listing will be made, 0.0 will get all, 999, for none). The information for 
components is requested for each unique component set, the three components of the fuselage were joined 
together because they were consecutively numbered (101, 102, 103 see files listing). Note also that the 
routine has automatically ignored interference shells (as will the wave drag) since they are not pertinent 
vehicle components . The viscous and wave drag programs also ignore jet flaps and slender bodies so the user 
can combine all the necessary components in the work file to do wave drag, viscous drag, and lifting surface 
analysis at one time. 
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BAY 

UOLUNE 

I 

117651.7 

2 

354885.9 

3 

87975.8 

4 

189873.8 

S 

41168.7 

s 

538.4 

7 

8589.4 

8 

2188.4 

s 

5849.7 

18 

13188.8 

li 

13638.1 

TOTAL 

747947.4 


iBSCXSSA SCALE 189.8 IN/1LOCK 
ORDINATE SCALE 1899.9 IN2/BLOCK 
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* 0 EXIT t i $F DRAG 8 2 EDIT 
? 

I 

SKIN FRICTION DRAG INPUT 
? ENT£R SAND GRAIN HEIGHT OCS(FT)) 

0*000022 

,® F0R AU * 9## fm nome > 

ENTER PRESS, TEnP,NACH(DI$P) FOR CASE i 

ENflR^RElsJf?^,nACH(DISP ) FOR CASE 2 
1455.8,483.0*1.2 

ENTER PRESS.TEHP,I1ACH(DISP) FOR CASE 3 
973.30*447.4.1.8 

ENTER PRESS, t£NP,NACH(DISP) FOR CASE 4 
? 


0 , 0.0 

FOR FUSELAGE 
7 

0 . 01,2 
INBD UING 
7 

0.01.1.2.2 
^OUTlD UING 

0*01,1.2*2 
^LOUER TIP FIN 

0.01,1.2.2 
^UPPER TIP FIN 

0.91,1.2.2 
LOWER UERT 

0.01.1.2.2 
VERT PYLON 

0.01.1.2,2 

CANARD 

7 

0 . 01 , 2 . 0,2 


ENTER! XTRANS/LENO, FLAT PLATE 1 t AXIS-0ODV 2 * 
ENTER! XTRANS/LENO, CK, FLAT PLATE 1 8 AIRFOIL 2 * 

ENTER! XTRANS/LENO, CK, FLAT PLATE 1 * AIRFOIL 2 t 

ENTER! XTRANS/LENO, CK* FLAT PLATE 1 8 AIRFOIL 2 8 

ENTER! XTRANS/LENO, CK, FLAT PLATE 1 8 AIRFOIL 2 8 

ENTER! XTRANS/LENO, CK, FLAT PLATE 1 8 AIRFOIL 2 8 

ENTER! XTRANS/LENO. CK, FLAT PUTS 1 8 AIRFOIL 2 8 

ENTER! XTRANS/LENO, FLAT PLATE 1 8 AXIS-IODY 2 8 
ENTER! XTRANS/LENO, CK, FLAT PLATE 1 8 AIRFOIL 2 8 
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* • EXIT t 1 $F DftM t 8 EDIT 
? 


After the viscous solution is concluded the user has the 
opportunity to edit the cross sectional data used in the 
calculations and/or add additional solutions to those already 
calculated. 

In this instance the user chose to end and the program proceeds 
into the next analysis routine selected (wave drag in this case) 
or returns to ANALY** if analysis is completed. 

Although each supersonic Mach number requested is displayed in the 
manner shown in the next two figures only M = 1.2 is shown here 
to add some compactness to this manual. Cross sectional area cuts 
as a function of roll angle are displayed first (ever fourth roll 
angle calculated is displayed ) and then a listing of drag as a function 
of roll angle, total drag and drag coefficient are printed. 






*Wlvtt Wave drag is always the last program to be executed and on completion 

nOKtt returns the user to ANALY**. The user may add other solutions, 

fore** eliminate an incorrect solution or in this instance exit and obtain 

displays of the force and drag results stored in the plot file. 

In the next figure the title sheet accompaning each call to FORCES is 
shown illustrating the solutions available, interpolated Mach numbers, 
t solution Mach numbers, reference data, and viscous solution input data. 

If a lifting solution is also stored control surface solutions, if any, 
are indicated, and whether the drag data is trimmed or untrimmed. 
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FORCt* 

ccUi 


ON 

oo . 


The keyword for FORCES is FORC**. The user selects one set of 
plots at a time and always returns to FORC**. 


The user first selected wave drag versus Mach number, then 
viscous drag, which exhausted the solutions stored in the 
plot file. 
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The user enters ANAIASIS and is told the previous 
solutions (Wave Drag) and the Mach numbers used 
( 0 . 6 , 1 . 2 ). 

The user elects to retain the same Mach numbers and 
asks for the lifting set up then goes into 
ANALYSIS program calculation (cj . 
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t CALCULATE LIFTING NATRIX t I Vff S i NOU 
t CALCULATE THICKNESS NATRIX tlVUtl NO<t 


» ENTER 3 DRAG CALC. ANGLES 
T 

1 I It 


t INITIAL CONRONENT INCIDENCE ANGLE t 


> INRUT WHEN CURSOR STORS UNDER INITIAL DEFLECTION t 


CONR. NO. CORR. NAME INITIAL DEFLECTION 


i 

2 . 

3 

4 

5 



The lifting set up is a set of numerical answers to 
questions posed by the computer and is organized to 
provide the user with a record of his input. 

The graphics cursor appears at the end of the set up 
to give the user a chance to review his input for. 
correctness and take a copy if he wishes. 
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NORflAL JOB COMPLETION 

COrWAND- purg«,ttp«l 
COWHAND- purg«,tap«3 

Cowhand- aliU g# Up«i l |i^ru #l d*fttm,rp«9M,M.oM II ,wfiU. 
NEU CYCLE CATALOG, CY-Ml 

COWHAND- c.UUg.Up*3* F UtfiU # id.it##M*rp-9M.cA-op.» # .x.fiU. 
NEU CYCLE CATALOG, CY-OOt 

COWHAND- c«ta log, Up»U , p»n« If 1 |* # t , rp«999 , 

NEU CYCLE CATALOG, CY-BBl 
COWHAND- logout 


Following job conpletion the user replaces* the old 
permanent geometry and plot files with the updated 

files by first purging the old files and then 
cataloging them again as new files. Some CDC 
facilities do not operate with this type of file 
management system and different methodologies for 
updating the Mass Storage files used in TAPE1 
and TAPE3 will be required. 

The user can submit a batch job to run the slender 
body and lifting surface solutions using the data 
listing on the following figure. 



100=PEDAK,T200,IO200,CM120000,P4. 

110=ACCT (DIVAN 0018422401*0111013) 

120 ss KEQOEST,TAPEl5 , *PF. 

130-REQUEST,TAPE6 , *Pf . 

140=ATTACH,TAPE4 , PANELFILE , ID=D0500. 
150==ATTACH,SLER,SLENLGAD,ID=D0500. 

160=MODE,0. 

170 = ATTACH , TAPE 2 , SC I FFGEOM , I D=D 0 500. 
180*REWIND,SLEN. 

190=SLEN. 

2 0 0=RETURN , SLEN , TAPE 4 , TAPE 2 . 

210=ATTACH # TAPE7 , PAN ELF I LE , ID=D0500 . 
220=ATTACH,TAPE14,PLOTFILE,ID=D0500. 

230 =ATTAC H , U DP , UDPLOAD , I D=D 0 500. 
240=REWIND,UDP. 

250=UDP • 

260=RETURW,TAPE9,TAPE10,TAP£11,TAPE12,TAPE13. 
270=CATALGG,TAPE6,UDPOUTPUT,ID=D0500,RP=20. 
280=CATALOG,TAPE15 , WAVEIN, ID=D0500 , RP=20 . 
290=CATALOG,TAPE14,PLOTFILE,ID=D0500,RP=20. 
300=EXIT. 

310=DISPOSE,TAPE6,PR=IAJ. 


, / 

The listing shown here was designed to execute the slender body solution 
then the lifting surface solution. If no slender body is present the slender 
body program automatically stops and goes directly to the lifting solution 
so this one control listing will perform all the necesary solutions. 

The output products of the program execution’s are, a detailed listing of 
program output (UDPOUTPUT on line 270), the wave drag due to lift pressure data 
file (WAVEIN) on line 280) and the updated plot file (PLOTFILE on line 290). 
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In this example the user is already into his session when he 
decides to investigate the solutions already present in the 
plot file, in this case for a single engine trainer concept. 


By enter ANALYSIS he finds that solutions exist for lifting 
surface, wave drag and that the isolated body solution was 
run with the lifting surface. By executing wave drag the user 
includes the transonic drag rise equations in his total drag solution. 


T v e user exits- from AFALYSIS and enters FORCES tc. see displays 
of stability derivatives and drag. The next figure shows the 
summary sheet which is always displayed at the start of the 
FORCES function. 
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FORCSS 

c«cl # cl* # cat 


The keyword for the FORCES function is 


FORC** 


The user requests displays of 


dc m/dc L 




All of which will be displayed versus Mach 
number, shown in the following figure. 
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The user now requests a display of induced drag versus lift for 
0.60 Mach number and the second suction variation stored in the 
Program (Cdi 2) 
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APPENDIX IV 


SYSTEM SUBROUTINE DESCRIPTIONS 
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ANALY: Set-up and executive to ANALYSIS function. Interprets user inputs, then 
schedules analysis programs or set-up program for sequential execution. 

AREA, AREAS ,AREAW: Similar routines designed to deteimined cross sectional cuts 

for selected Mach angles and roll angles. AREA and AREAS are preset to Mach - 1.0, 
ROLL - 0.0 degrees. 

COD, CODIM: Second order curve fit routines which appear in different segments of 

the program structure. 

CPLOT: Provides a simple unlabeled ordinate versus abscissa plot for specified 

number of input points 

CPFR: Set up for spanwise and chordwise plots. Sums combinations of user 

input parameters (a, 3, P, Q , R, etc) and display chordwise or spanwise loads for 
those conditions. Also displays twist, camber and thickness distributions. 

CPWING: Summing routine called by CPPR to calculate chordwise increments. 

DECRD: Card image numerical processing routine. Interprets card image values 

for correctness and fills input array leaving blank locations unchanged. 

DIGIT: Digitizing routine for cross-sectional input. Builds new components with 

sequential sections or re -digitizes old components with arbitrarily ordered 
sections . 

DISP: Orthographic projection of two or more components in a file. 

DISPS: Orthographic projection of single component. Also displays pertinent data 

for component and allows adjustment of longitudinal display lines. 

DQWS: Calculates total wave drag for input configuration, prints out result, and 

stores total drag coefficient. 

EDITN: Body component visual editor. Displays component to be edited and 
uses graphics cursor to edit component. 

EDITP: Surface component visual editor. Functions as spanwise editor and is used 

for paneling components for lifting surface solutions. Principle method of entry 
is graphics cursor. 

EDPUT: Executive routine for LIST function. Interprets alphanumerics so LIST can 

perform the indicated function. 

EMLORD: Calculates Eminton-Lord isolated body optimum area distribution at 
selected output stations. 
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FIT: Checks input curves for breaks and splines or least squares individual 

segments . 

FITWT: Least squares curve fitting routine designed with fixed endpoints and 

special polynomial leading edge fit for airfoils with leading edge radius. 

FORCE: Set-up for total coefficient and drag coefficient display. 

FOREX: Executive for FORCE, interprets user inputs into commands and numerical 

values . 

GETPNT : Modified Tektronix routine for inputing a point off the digitizing 
tablet. 

GEOEX. Principle executive interpreter. Reduces alphanumeric and numeric inputs 
into control numerics for MAIN program execution. ' 

GREEK: Produces graphical constructions of A, r and A for TYPE 3 and 4 title block 

under DISPS . 

GRID: Provides five standard grids and numbering used under subroutine PLOTIT 

HSHLDR: Linear least square equation solution routine. 


ICON: Simple geometry menu sub -program. Request numerical inputs, then routes 

program to routines designed to intemret request. 

IDATA: Viscous drag input -menu sub-program. Request data for cases and completes 

geometric definition for viscous solution. 

INFO: Calculates and displays component data out of DISPS. 

INSTR1 : Menu display for EDIT functions. 

INSTR2 : Menu display for INTERFERENCE and THREE-VIEW functions 

INTERF: Graphical input of interference shells and jet flaps using graphics cursors. 

LABEL, LABEL7: Label subroutines for plots from LOAD and FORCE functions 

respectively. 

LIMIT, LIMITS, LIMITW : Calculates component limits needed for subroutines AREA, 

AREAS, and AREAW respectively. 

LIST: Cross section and keyboard editing routine. Display true view of cross 

sections and/or point listing. 
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LSET: Converts surface components into panel description. Ques user for flap 

data, and initial data concerning lifting surface solutions. 

MATRIX: Sets up boundary conditions and solves EMLORD coefficients. 

MAXIT: Determines parametric minimums and maximums under subroutine PLOTIT. 

MJLPNT: Combines with GETPNT for multiple point digitizing from tablet. 

ORIG:. Places a graphical square at specified points on the screen. 

PLEXC: Executive interpreter for CPPR subroutine. 

PLOTIT: Graphical plot routine for LOADS and FORCES function. Will plot data on 

five grids. Axis names have been pre-set in the routine. 

PLOTOT: Simple display of planar components for EDIT function. 

SAVEB: Interprets and curve fits digitized data from subroutine DIGIT. Prompts 

user for necessary information. 

SCALIT: Selects scaling values used is subroutine PLOTIT. 

SDQ: Calculates wave drag for single roll angle. 

SECT: Utility routine for section deleting and interpolation; all component 

types. 

SHEET: Displays the title sheet at the initiation of FORCES function. 

SHI FT, SHIFTS, SHI FTW: Reorganizes points in cross section cuts made by AREA, 

AREAS, and AREAW for proper integraion. 

SLEND: Routine which combines .body components to form one slender body. 

SPLINE: Third order open ended spline routine, used as principle routine for 

section fitting, curve integrations and differentiation. 

SYM4,SYM: Symbol making routines used in various sections of program. 

TABARM: Turns terminal control from the keyboard to the tablet. 

TABINT: Initializes tablet variables for tablet input. 

TABOFF: Turns off digital tablet control. 

TLU: Linear Interpolation routine used in viscous flow solutions. 
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TOTAL: Interprets output from analysis programs, adds empirical data where 

necessary, and stores the results. 

TRANS. Inputs component translations from keyboard or graphics cursors, then 
translates components. " ? 

TRAP: Calculates simple wing geometry and sets up data to be converted to 

component geometry formats . 

VDRAG: Viscous drag set-up routine. 

VIEW3: Three view projection routine used for three views and interference 

shells . ' 

WCAD: Card image component interpreter. 

WDRAG : Main set-up routine for wave drag calculation. 

WET: Calculates wetted area and volume for components in viscous drag routines. 

WSVP: Viscous drag solution routine. Uses data from subroutine. 

IDATA. Prints and stores answers. 

XACALW: Computes minimums 1 and maximums for combined components in wave drag. 

XYZ,XYZR: Routines to re-organize component description from one form to another 

similar more descriptive form. XYZ converts card data to the component 
descriptions, XYZR changes geometry for wave drag calculations. 
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